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Notation 

Cp  -  Pressure  coefficient  defined  in  eq.  (3) 

a  -  Yaw  angle  relative  to  laboratory  space 

<j>  -  Pitch  angle  relative  to  laboratory  space 

-  Total  pressure 

P  -  Static  Pressure 
s 

arp  ”  Yaw  ang^e  rolative  to  probe  axis 
<j>  rp  -  Pitch  angle  relative  to  probe  axis 
Pp  -  Pressure  indicated  by  the  probe 
dp  -  Yaw  setting  of  the  probe 
4> p  -  Pitch  setting  of  the  probe 


Subscripts 


I 

II 

III 

IV 
A 
B 

Lo 

rp 

Up 


I 

I  relating  to  probe  positions  shown 
in  figure  3 

relating  to  A  type  probe 
relating  to  B  type  probe 
lowest  value  in  an  array 
relative  to  the  probe  axis 
highest  value  in  an  array 
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Notation  of  PENPTS 


DO  -  denominator  In  the  subsequent  line  of  the  program 

HT  -  penetration  height  (Cp  of  a  calibration  surface) 

IR  -  a  flag  with  a  value  of  either  1  or  2  to  identify  first 

or  second  penetration  point.  This  flag  also  controls 
return  of  values  to  calling  program  or  subroutine 

IB  -  a  flag  indicating  the  location  of  the  last  rectangle 
in  the  band  checked  for  a  penetration  point.  It  pre¬ 
vents  repetition  of  the  scanning  in  the  103  Do  loop 

NX  -  number  of  X  values  in  the  calibration  surface  grid  (rvalues  of  the 
calibration  surface) 

NY  -  number  of  Y  values  in  the  calibration  surface  grid  (.  values  of  the 
calibration  surface) 

X(I)  -  X  values  in  the  grid 

XM  -  value  of  X  at  point  B  (fig.  5) 

XR( I)  -  two  values  of  X  returned  by  the  subroutine  as  results 

XS  -  currently  calculated  value  of  X  which  is  the  present 
result  and  is  later  stored  in  XR(I) 

Y ( J )  -  Y  values  in  the  calibration  grid 

YC  -  intermediate  value  used  in  the  subsequent  line  of  the  program 

YG  -  the  Y  position  of  the  penetrating  line* 

Z(I.J)-  Values  of  the  Z  =  Cp  calibration  surface  above  grid  points 

ZX  -  value  of  Z  either  at  point  A  or  point  C  (see  fins.  5, 

6  and  7) 

ZM  -  value  of  Z  at  point  B 


*  It  should  be  pointea  out  here  that  the  penetrating  line 
is  parallel  to  the  X  axis  at  a  height  HT  =  Cp  above  the 
X,Y  plane  and  at  a  distance  YG  from  the  X  axis  on  the  X,Y  plane 
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Notation  of  INTSCS 


AA  -  constant  defined  in  eq.  8 

AB  -  constant  defined  in  eq.  9 

ARIN  -  ordinate  at  which  scanning  starts 

ARLO  -  lowest  ordinate  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  all  calibration 
surfaces  used) 

ARM1  -  ordinate  measured  relative  to  probe  1  (the  probe 

belonging  to  the  first  curve  Cp  =  const.) 

ARM2  -  ordinate  measured  relative  to  probe  2  (the  probe 

belonging  to  the  second  curve  Cp  =  const) 

ARP1  -  setting  ordinate  of  probe  1 

ARP2  -  setting  ordinate  of  probe  2 

ARRES  (1)  -  abscissa  value  returned  to  INTSCS  from  PENPTS 

ARUP  -  highest  ordinate  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  all  calibration 
surfaces  used) 

AR2R  -  ordinate  relative  to  laboratory  space 

AR2RJ  -  ordinate  relative  to  laboratory  space  of  previous 

scan 

BA  -  constant  defined  in  eq.  8 

BB  -  constant  defined  in  eq.  9 

DAR  -  ordinate  step  for  scanning 

HT1  -  penetration  height  of  calibration  surface  of  probe 

1  (Cp  of  probe  1) 

HT2  -  penetration  height  of  calibration  surface  of  probe 

2  (Cp  of  probe  2) 

-  a  flag  which  is  equal  to  2  when  the  calculation  is 
carried  out  as  shown  in  fig.  11,  otherwise  it  is 
equal  to  1.  (This  information  is  required  in  the 
main  program  and  is  returned  to  it) 


IEPS 


ISL  -  a  flag  which  is  equal  to  10  if  one  of  the  left 

hand  side  results  returned  from  PENPTS  is  1000.0 
(i.e.,  is  not  a  penetration  point).  Otherwise 
it  is  equal  to  1 

ISR  -  a  flag  which  is  equal  to  10  if  one  of  the  right 

hand  side  results  returned  from  PENPTS  is  1000.0 
(i.e.,  is  not  a  penetration  point).  Otherwise 
it  is  equal  to  1 

K  -  a  flag  indicating  whether  the  first  or  the  second 

intersection  is  evaluated 


Nl  -  number  of  abscissas  in  calibration  surface  grid 

(fig.  4) 

N2  -  number  of  Ordinates  in  calibration  surface  grid 

(fig.  4) 

RES1(K)  -  evaluated  and  returned  abscissa  of  intersection 

RES2(K)  -  evaluated  and  returned  ordinate  of  intersection 

R1  -  radius  in  fig.  11  for  probe  1 

R2  -  radius  in  fig.  11  for  probe  2 

RAIL  -  abscissa  of  point  AIL  in  fig.  9  or  10 

RAIR  -  abscissa  of  point  AIR  in  fig.  9  or  10 

RAJL  -  abscissa  of  point  AJL  in  fig.  9  or  10 

RAJR  -  abscissa  of  point  AJR  in  fig.  9  or  10 

RBIL  -  abscissa  of  point  BIL  in  fig.  9  or  10 

RBIR  -  abscissa  of  point  BIR  in  fig.  9  or  10 

RBJL  -  abscissa  of  point  BJL  in  fig.  9  or  10 

RBJR  -  abscissa  of  point  BJR  in  fig.  9  or  10 

XI ( I)  -  abscissas  of  calibration  surface  of  probe  1  or 

Its  transformation  as  determined  by  calling 
statement  in  main  program 

X2( I )  -  abscissas  of  calibration  surface  of  probe  2  or 

its  transformation  as  determined  by  calling 
statement  in  main  program 
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X3(I)  -  abscissas  of  calibration  surface  of  probe  1  or  its 

transformation  as  determined  by  calling  statement 
in  main  program 

X4(I)  -  abscissas  of  calibration  surface  of  probe  2  or  its 

transformation  as  determined  by  calling  statement 
in  main  program 

Y1(I)  -  ordinates  matching  XI (I) 

Y2(I)  -  ordinates  matching  X2(I) 

Y3(I)  -  ordinates  matching  X3(I) 

Y4(I)  -  ordinates  matching  X4(I) 

Z1(I,J)  -  values  of  Cprp  for  XI (I)  and  Y1(I) 

Z2(I,J)  -  values  of  Cppp  for  X2(I)  and  Y2(I) 

Z3(I,J)  -  values  of  Cprp  for  X3(I)  and  Y3(I) 

Z4(I,J)  -  values  of  Cprp  for  X4(I)  and  Y4(I) 
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ALFCN 


Notation  of  Main  Program  VDR 


-  variables  in  library  routine  IXCLOK  used  to  evaluate 
computation  times 

-  values  of  a  returned  to  VDR  from  INTSCS 

-  a  input  when  experiment  simulation  is  carried  out 


rp  III 

new  value  of  'for  next  iteration 


ALFDN  -  new  value  of  arpjy  for  next  iteration 


-  mach  number 


EPSPS  - 


VpIII 


relative  difference  between  present  static  pressure 
and  its  value  in  the  previous  iteration 


EPSPSG  -  convergence  criterion  (on  the  static  pressure) 

FALF  -  final  a  value  (at  convergence) 

FPHI  -  final  4  value  (at  convergence) 

IA50  -  a  flag  in  IXCLOK 

IA60  -  a  flag  in  IXCLOK 

I  COPS  -  number  of  PS  corrections  carried  out  to  ensure  the  Cp's 
are  in  calibration  range 

IEPS  -  a  flag  governing  the  convergence  criterion  returned 
from  INTSCS 


IEPS1  -  a  flag  governing  the  convergence  criterion  for  probes 
I  and  II 

IEPS2  -  a  flag  governing  the  convergence  criterion  for  probes 
I  and  III 

I  IT  -  number  of  iterations  on  ps 

ISCAN  -  number  of  static  pressure  scans  from  initial  static 
pressure  guess  upwards 

IXCLOK  -  system  subroutine  for  computation  time  evaluation 

NOCOPS  -  a  flag  which  when  equal  to  1  causes  static  pressure 
corrections  to  be  carried  out  such  that  Cp's  are 
always  inside  calibration  range,  and  when  equal  to  2 
cause  these  corrections  to  be  skipped 

NOSIM  -  a  flag  which  when  equal  to  1  causes  the  program  to 
simulate  velocity  measurement  experiments  and  when 
equal  to  2  causes  the  program  to  reduce  measured  data 

NX  -  number  of  calibration  surface  matrix  abscissa's, 
a  values  (must  be  identical  for  all  calibration 
surface  matrices  used) 

NY  -  number  of  calibration  surface  matrix  ordinates, 

<j>  values  (must  be  identical  for  all  calibration  surface 
matrices  used) 

>DYN  -  dynamic  pressure 
I  -  value  of  <j>  returned  to  VDR  from  INTSCS 


>HIS  -  %u 

‘HI I  -  $  input  when  experiment  simulation  is  carried  out 
>HI1  -  ♦rpl 

>HI2  '  VpII 

>HI3  -  *rpl 

’n14  ‘  Vpin 

>PA  -  signal  of  probe  in  position  I 

>PB  -  signal  of  probe  in  position  II 


PPC  -  signal  of  probe  in  position  III 

PPC  -  signal  of  probe  in  position  IV 

PS  -  static  pressure 

PS  IN  -  initial  guess  of  static  pressure 

PSN  -  new  value  of  static  pressure  for  next  iteration 

PSC  -  corrected  static  pressure 

PST  -  memorized  corrected  static  pressure 

PT  -  total  pressure  measured 

RELXPS  -  relaxation  factor  for  convergence  on  the  static 
pressure 

VIRTIM  -  virtual  computation  time 

XA  -  abscissas  of  calibration  surface  matrix  of  probe  A 

XAX  -  transformed  abscissas  of  calibration  surface  matrix 

of  probe  A 

XB  -  abscissas  of  calibration  surface  matrix  of  probe  B 

XBX  -  transformed  abscissas  of  calibration  surface  matrix 

of  probe  B 

XPC  -  yaw  setting  of  position  III 

XPD  -  yaw  setting  of  position  IV 

XRIN  -  a  value  at  which  scanning  starts 

XRLU  -  lowest  a  value  of  calibration  surface  (as  a  condition 

it  must  be  identical  for  probes  A  and  B) 

XRUP  -  highest  a  value  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  probes  A  and  B) 

YA  -  ordinates  of  calibration  surface  matrix  of  probe  A 

YAX  -  transformed  ordinates  of  calibration  surface  matrix 

of  probe  A 

YB  -  ordinates  of  calibration  surface  matrix  of  probe  B 

YBX  -  transformed  ordinates  of  calibration  surface  matrix 

of  probe  B 
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YPA 


-  pitch  setting  of  probe  A 

-  pitch  setting  of  probe  B 

-  <j>  value  at  which  scanning  starts 


YPB 
YRIN 

YRLO  -  lowest  $  value  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  probes  A  and  B) 

YRUP  -  highest  4  value  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  probes  A  and  B) 

ZA  -  CprpA  va1ues  of  Probe  A 

ZAX  -  transformed  CprpA  values  of  probe  A 

ZAMAX  -  maximum  of  ZA  array 

ZAMIN  -  minimum  of  ZA  array 

ZB  -  CprpB  va1ues  of  Probe  8 

ZBX  -  transformed  CprpB  values  of  probe  B 

ZB MAX  -  maximum  of  ZB  array 

ZBMIN  -  minimum  of  ZB  array 


9 


Introduction 


Experimental  knowledge  of  the  flow  field  generated 
by  rotating  turbo  impellers  is  of  prime  importance  in  the 
research  and  development  of  turbomachinery.  It  is  essential 
for  the  refinement  of  design  methods,  for  the  development 
of  new  flow  models  which  include  secondary  flow  and  tip 
clearance  effects,  and  particularly  for  the  verification 
of  new  computer  codes  developed  to  calculate  the  flow  through 
rotating  blade  rows. 

In  recent  years  laser  velocimeter  techniques  have 
been  applied  successfully  to  measure  the  flow  both  inside 
and  downstream  of  rotors.  (Ref.  1  for  example).  It  has 
become  clear  however,  that  the  laser  techniques  are  only 
reliable  in  the  hands  of  experienced  i nvesti gators .  A 
window  which  remains  clean  is  essential,  and  seeding  is 
usually  required.  Laser  techniques  do  not  measure  the 
pressure  field  and  usually  can  only  measure  two  components 
of  the  velocity  unless  the  axis  of  the  laser  is  tilted. 
Difficulty  is  also  encountered  when  measuring  close  to 
walls.  Hence  there  are  reasons  to  consider  alternative 
techniques,  particularly  if  they  are  simpler  to  apply  routinely 
in  stationary  turbomachinery  passages.  Furthermore,  the 
achievement  of  redundancy  in  measuring  the  flow  field 
behind  the  impeller  is  itself  a  worthwhile  goal.  The  present 
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work  deals  with  the  application  of  a  particular  system 
of  small  high  response  pressures  probes  at  the  exit  of 
an  Impeller. 

Measurements  behind  an  impeller,  in  the  stationary 
bladeless  gap,  are  simpler  to  make  than  measurements  within 
the  rotating  passages.  Transducer  probes  can  be  installed 
through  the  stationary  machine  casing  and  the  data  transmitted 
without  resort  to  slip  rings  or  rotary  transformers.  The 
sensor  is  not  subjected  to  the  centrifugal  field  or  to 
the  vibration  of  the  rotor.  However,  the  flow  to  be  measured 
is  then  fluctuating  at  blade-passing  frequency  and  any 
system  of  sensors  must  be  calibrated  for  a  wide  range  of 
possible  Mach  number,  pitch  angle,  yaw  angle  and  pressure 
variation  -  and  yet  must  be  capable  of  the  necessary  fre¬ 
quency  response. 

In  Ref.  2,  a  method  was  described  for  using  two  semi¬ 
conductor  pressure  probes  together  with  the  technique  of 
synchronized  sampling,  to  obtain  the  distribution  of  the 
velocity  vector  downstream  of  a  rotor.  The  geometries 
of  the  two  probes,  designated  Type  A  and  Type  B,  and  their 
installation  in  the  compressor  annulus  are  shown  in  Fig 
1.  It  was  argued  that,  in  principle,  by  rotating  the  probes 
in  yaw  about  their  tips  and  controlling  the  sampling  of 
the  data  from  each  probe  to  be  at  the  same  position  in 
the  rotor  frame,  the  system  of  two  separate  probes  could 
be  used  to  acquire  data  at  a  point  in  the  periodic  flow 
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from  the  rotor,  corresponding  to  data  normally  obtained 
from  the  multiple  sensors  of  four-  or  five  hole  pneumatic 
probes  when  measuring  velocity  in  steady  flows.  The  tech¬ 
nique  promised  the  use  of  probes  having  the  simplest  geo¬ 
metry  and  thus  avoided  the  large  size,  expense  and  unreli¬ 
ability  of  multiple  sensor  probes  which  incorporate  multipl 
semi-conductor  transducers  (Ref.  3).  Because  of  the  simple 
sensor  tip  geometry  (that  of  a  cylinder  at  incidence  to 
the  flow)  the  unsteady  response  was  likely  to  be  as  good 
as  could  be  expected  of  any  single  physical  sensor. 

The  two-probe  technique  is  strictly  applicable  only 
to  periodic  flows.  However,  data  obtained  on  successive 
rotations  of  the  rotor  can  be  averaged  to  eliminate  fluctua 
tions  which  are  not  periodic.  This  was  shown  to  be  effect¬ 
ive  in  tests  reported  in  Ref.  2  in  which  a  single  Type 
A  probe  was  used  to  establish  the  peripheral  bl ade- to-bl ade 
distribution  of  flow  yaw  angle. 

In  order  to  obtain  velocity  from  the  pressure  measure¬ 
ments  which  can  be  obtained  from  the  two  probes,  the  steady 
response  characteristics  must  first  be  established  in  cali¬ 
bration  tests  carried  out  in  a  known,  controlled,  uniform 
flow.  Second,  a  method  of  applying  the  calibration  to 
measurements  made  in  an  unknown  flow  must  be  devised. 

In  the  present  method  two  different  approaches  have  been 
followed.  In  Ref.  4,  a  technique  is  given  for  representing 
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and  applying  the  probe  calibration  analytically.  When  first 
applied  the  method  gave  surprisingly  good  accuracy  (1  -  2$)  since 
the  method  required  that  the  probes  had  characteristics 
which  could  be  well  represented  analytically.  This  in 
turn  required  that  the  probe  tips  be  geometrically  precise, 
a  feature  which  was  not  present  in  the  first  generation  of 
probes . 

A  second,  more  general  approach  is  reported  here,  wherein 
the  calibration  of  each  probe  is  represented  by  a  two-dimen¬ 
sional  array  of  pressure  coefficients.  The  application  of 
the  calibration  given  in  this  form,  in  an  unknown  flow  required 
the  development  of  special  numerical  procedures.  The  purpose 
of  the  present  report  is  to  document  the  analysis  and  the 
Fortran  program  developed  to  apply  the  method. 

In  its  present  form,  the  method  does  not  require  that 
the  calibration  "Surfaces"  be  symmetrical  about  any  axis  or 
be  expressed  in  analytical  form,  but  does  require  that  the 
pressure  coefficient  be  independent  of  Mach  number.  The  lat¬ 
ter  restriction  could  undoubtedly  be  removed  by  introducing 
additional  iterative  steps.  Further,  in  the  present  method 
only  five  measurements  have  to  be  taken  to  determine  uniquely 
a  velocity  vector  at  a  point.  Throughout  the  report,  the 
Fortran  program  notation  has  been  used  to  describe  the  physics 
and  equations  involved  in  the  solution. 


Mathematical  model 


Assume  the  A  and  B  type  probes  of  ref.  [1]  (see  also  fig.  1) 
to  be  immersed  in  a  three-dimensional  steady  flow  field.*  The 
pressure  response  of  each  of  these  probes  in  given  gas  is  func¬ 
tionally  described  by  four  variables  as: 

Pp  *  %  (a.*  ,  PT.  P$)  (1) 

If  a  pressure  coefficient  is  defined  as 

p  -  p 

c%  -,vni7  (2) 

The  calibration  surface  of  each  probe  is  given  in  the  general 
case  in  form  of  a  matrix  of  values  of  Cp,  where 

Cprp  =  Cprp(arp’^rp)  (3) 

The  pressure  coefficient  defined  in  this  way  has  only  a  second 
order  dependence  on  Mach  and  Reynolds  numbers  in  the  range  of 


*  For  our  purpose,  using  the  synchronized  sampling,  the  flow 
field  behind  the  impeller  is  steady,  although  the  probes  re 
quire  a  high  speed  response  because  of  fluctuations. 
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0  <Ma<  0.7  in  turbulent  flows  (Ref.  2),  so  that,  to  first 
order  their  influence  is  neglected  in  writing  eq .  (3). 

If  the  type  A  probe  is  rotated  about  its  axis  into  three 
different  positions  (i  =1,  III,  IV)  readings  are  taken,  and 
the  type  B  probe  is  fixed  in  position  II  and  a  single  reading 


is  taken,  the  following  four  equations  can  be  written, 


CPAi  ‘ 

'PB 1 1  '  Ps  ’  CpBI 


it  should  be  pointed  out  here,  to  avoid  misunderstanding  that  a 
and  4>  are  defined  in  a  coordinate  system  relative  to  the 
machine  axis  and  not  relative  to  the  probe  axis.  In  the  set 
of  four  equations  (4)  there  are  four  unknown  quantities,  namely: 
o,  <j>,  Pj  and  P$  .  These  are  the  quantities  to  be  evaluated  using 
the  measured  data.  Together  with  the  stagnation  temperature  they 
define  the  flow  field  uniquely.  The  four  equations,  resulting 
from  the  four  measurements,  should  be  sufficient  to  determine 
the  four  unknown  quantities. 

But  the  problem  is  complicated  by  three  facts: 

1)  The  calibration  surfaces  are  not  generally  known  in 
analytical  form. 

2)  The  calibration  surfaces  are  double  valued  in  a  and  $ 
i.e.,  for  a  given  Cp  and  a  there  exist  two  41 
values,  or  for  a  given  Cp  and  41  there  are  two  a 
values  which  satisfy  eq.  3. 


16 


* 

/ 


3)  As  a  result  of  the  measurement  the  value  of  P  rather  than 

P 

of  C  is  determined. 

P 

Since  the  calibration  surfaces  are  not  given  in  a  simple  analytical 

form  the  solution  has  to  be  numerical.  An  iterative  procedure  is 

required  because  P_  and  not  C  values  are  measured.  First  P  and 

P  P  s 

P.J.  have  to  be  guessed  to  yield  Cp  values,  knowing  the  four  measured 
values  of  Pp,  and  using  eq.  4.  The  guess  is  then  iteratively 
corrected  to  converge  on  the  solution.  However,  convergence  of  the 
iterative  procedure  is  complicated  because  of  the  double  valued 
nature  of  the  calibration  surfaces. 

This  method  of  iteration  shown  in  Fig.  2  was  attempted  initially 
for  the  evaluation  of  a  measured  point.  In  practice  the  initially 
suggested  procedure  converged  in  some  cases  and  diverged  in  others, 
depending  on  the  values,  and  signs  of  a  and  <f>.  This  was  not 
surprising  as  convergence  on  two  variables  is  not  likely  to  be  a 
simple  matter.  However,  in  order  that  the  measurement  technique 
be  useful,  convergence  on  the  correct  solution  for  a  general  set  of 
measurements,  is  absolutely  necessary.  In  practice,  this  can 
certainly  be  achieved  if  one  of  the  two  iteration  variables  is 
obtained  by  measurement.  Since  the  static  pressure  is  a  difficult 
quantity  to  measure  even  in  a  steady  flow  field,  only  the  stagna¬ 
tion  pressure  measurement  need  be  considered.  It  is  possible  that 
the  time-varying  stagnation  pressure  could  be  measured  with  a  suit¬ 
ably  designed  Kiel  probe.  Data  would  then  be  taken  by  synchronized 
sampling  from  the  fixed  Kiel  probe,  from  the  type  A  probe 
rotated  into  two  positions,  and  from  the 
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g.  2:  Unstable  iterative  solution  using  measurement 
of  A  type  and  B  type  probes  only. 


fixed  type  B  probe.  The  method  of  solution  is  then  as  shown 
in  Fig.  3.  The  method  shown  in  Fig.  3  proved  to  converge  under  all 
conditions.  It  is  described  in  detail  in  the  following  pages. 
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Stable  iterative  solution  using  measurements 
of  A  type  and  B  type  probes  as  well  as 
measurements  of  a  Kiel  probe. 


ffl 


Evaluation  of  a  and  it  (method) 

The  best  way  to  understand  the  evaluation  of  yaw  and  pitch 
is  to  look  at  it  from  a  topographical  point  of  view.  Each  of 
the  type  A  and  type  B  probes  has  a  unique  calibration  surface 
Cprp  =  cPrp^arp’  Vp^  *  where  arp  dnd  $rp  are  measured  relative 
to  the  axis  of  symmetry  through  the  sensor  at  the  probe  tip. 

The  calibration  surface  in  this  representation  is  invariant  to 
yaw  and  pitch  of  the  probe  axis  relative  to  the  laboratory 
space.  The  same  calibration  surface  if  represented  in  the  form 
Cp  =  C  p  ( oc ,  <t> )  can  be  written  as 

Cp  =  C  p  [  ( a  p  +  arp),  Up  +  *rp)]  (5) 

where  ap  and  <f>p  are  the  probe  tip  axis  angular  settings 
relative  to  the  laboratory  space.  It  is  clear  from  eq .  (5) 
that  Cp  can  be  derived  from  Cprp  by  a  constant  translation: 

“p  *  'fp  on  t*ie  a  >  $  plane.  As  each  Cprp  can  be  viewed 
as  a  hill  with  its  peak  at  arp  =  0  and  <f>rp  =  0.  The  Cp 
surfaces  are  the  same  hills  with  their  peaks  translated  to 
cx  =  ap  and  $  =  <j>  p. 

In  the  present  method  probe  A  is  used  in  three  different 
angular  settings  namely: 

I)  ap  =  0  *p  =  0 
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m)  ap  =  aplll  *p 

IV)  ap  =  apl V  *p 


0 

0 


and  the  B  probe  is  used  in  a  fixed  position,  namely: 


For  this  case  the  topography  of  the  calibration  surfaces  will 
appear  as  four  hills.  The  three  with  their  peaks  at  points 
(0,0);  ( +ap i i i »  0)  and  (-apjV,  0)  are  the  translated  hills  CprpA 
and  the  fourth  with  its  peak  at  (0,  4>pjj)  is  the  translated 
hill  CPrpB*  Their  contours  of  constant  Cp  are  then  as  shown 
in  figure  4. 

Assume  now  that  a  velocity  vector  with  yaw  a  and  pitch 
<p  is  to  be  measured.  These  values  of  a  and  $  will  be 

sensed  uniquely  by  the  probes  at  their  four  angular  settings.  Were 
equation  (4)  single  valued,  the  values  of  a  and  <»  could  be 
uniquely  evaluated,  as  the  single  intersection  point  between  the 
projections  of  appropriate  lines  of  constant  Cp  on  the  a,$  plane*. 
In  the  present  case  of  the  double  valued  functions,  the  lines  of 
constant  Cp  are  closed  curves  and  more  than  a  single  intersection 
point  do  exist. 


*  Cp  has  generally  a  different  value  for  each  probe  in  each 
of  its  settings.  Thus  the  solution  involves  solving  for  the 
intersection  points  between  projections  of  contours  of  specified 
(but  different)  Cp  values  on  the  different  hills. 
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Projections  of  Cp  =  const,  lines  of  the  four 

calibration  surfaces  (the  center  hill  is  for 
probe  A  at  ap  =  0,  <{>p  =  0,  the  top  hill  is 

for  probe  B  at  op  =  0,  4>p  =  $pII,  the  right 

hill  is  for  probe  A  at  ap  =  +apIII,  <f>p  ■  0 

and  the  left  hill  is  for  probe  A  at  ap  =  apiv’ 

4>p  s  °. 


The  situation  is  shown  in  Fig.  5  for  an  example  of  a  number  of 
such  intersection  points.  The  correct  intersection  point,  however, 
is  uniquely  identified  as  the  only  point  through  which  all  four 
Cp  *  const,  curves  pass. 

From  an  examination  of  Fig.  5  it  is  clear  that  the  inter¬ 
section  points  of  three  closed  Cp  curves  projections  only  are 
sufficient  to  identify  a  and  $  uniquely.  However  one  of 
these  three  curves  must  be  that  belonging  to  the  type  B  probe. 

The  details  of  the  numerical  procedure  used  to  obtain  the 
correct  intersection  point  in  the  evaluation  of  a  and  $  , 

are  given  in  the  following  paragraph. 
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Evaluation  of  the  intersection  point  coordinates  (numerical 
Procedure ) 


The  calibration  surfaces  Cprp  are  represented  in  form 

of  a  linear  string  of  values  ordered  into  an  array,  as  shown 

schematically  in  figure  6.  The  numbers  in  Fig.  6  indicate  the 

position  of  a  Cprp  value  in  the  string.  The  strinq  starts 

with  the  value  of  Cprp  at  a  point  (-arpLo;  -i>rpLo) 

and  ends  with  the  value  of  Cp  „  at  a  point  (+«  . 

rrp  K  rpUp  rpllp 

This  sequence  must  be  kept  and  can  not  be  changed. 

Each  of  the  closed  Cp  =  const,  curves  projections  in  figs. 

2  and  3  is  the  projection  of  the  line  of  intersection  between  a 
plane  parallel  to  the  a ,  <p  plane  at  a  height  equal  to  Cp  and 
the  calibration  surface.  These  closed  curves  can  be  determined  as 
the  locus  of  the  projections  of  all  the  penetration  points  of  ar¬ 
bitrary  lines  parallel  to  the  a ,  <t>  plane  at  a  height  Cp  and  the  call 
bration  surface.  Such  penetration  points  are  calculated  in  sub¬ 
routine  "PENPTS",  Fig.  10. 

PENPTS  calculates  the  first  two  penetration  points  of  a 
surface  Z  (in  the  present  case  Z  =  Cp)  by  a  straight  line 
piercing  that  surface.  If  the  surface  is  double  valued  these 
two  points  are  the  only  roots.  The  surface  is  given  as  a  table 
of  numbers  on  a  cartesian  basis  Z  =  Z(X,Y)  (or  in  the  present 
case  Cp  =  Cp( o ,  <f) . 

The  subroutine  has  the  following  limitations: 

1)  No  roots  can  be  found  on  the  lower  Y  =  const,  boundary 
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2)  The  roots  can  be  found  only  along  a  piercing  line 


3) 


4) 

5) 

6) 

7) 


which  is  parallel  to  the  X  axis  (but  at  any  heiqht 
above  the  X,Y  plane) 

The  table  Z  =  Z ( X  ,  Y )  de  fining  the  surface  must  be 
based  on  a  grid  comprising  lines  X  =  const,  and 
Y  =  const.  The  spacing  of  these  lines  must  not  be 
equal.  In  other  words  the  surface  is  defined  by  a 
rectangular  grid  in  the  X,Y  plane,  from  X  .j  to 


X„,.„  and  from  Y  • „  to 
max  mi n 


max 


Only  a  single  root  can  be  evaluated  in  a  surface 
element  located  above  a  defining  rectangle  on  the 
X,Y  pi ane  . 

Not  more  than  the  first  two  roots  will  be  evaluated  for 
any  piercing  line. 

The  surface  must  be  monotonic  over  each  rectangle 
(this  is  a  result  of  limitation  4). 

All  X  and  Y  arguments  must  be  given  in  increasing  order. 


These  limitations  do  not  restrict  the  application  of  PENPTS 
in  the  present  problem  as  long  as  the  calibration  surfaces  are 
smooth  within  the  element  located  above  a  grid  rectangle. 

However,  the  elemental  grid  rectangles  can  be  reduced  a  rbi  - 
trarily  in  size.  If  a  calibration  surface  is  more  than  double 
valued  PENPTS  will  fail.  However  in  this  case  the  probe  yielding 
such  a  calibration  surface  can  not  be  considered  a  useful  instrument 
unless  used  only  in  parts  of  the  domain  where  it  is  double  valued. 
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Z  Is  defined  from  Xm.  to  Xma„  and  from  Y.„  to  Y„„„. 

min  max  ml n  max 

PENPTS  is  given  the  following  initial  information:  X  values  and 
Y  values  defining  grid  points,  corresponding  Z  values,  the  Y 
location  (Y6)  and  the  height  above  the  X, Y  plane  ( HT )  of  the 
piercing  line.  With  this  information  PENPTS  searches  for  the  band 
of  rectangles  which  includes  YG  or  of  which  YG  is  the  lower  bound¬ 
ary  (see  fig.  7)  and  then  scans  this  band  from  left  to  right  in 
search  for  penetration  points.  The  scanning  is  based  on  the  geo¬ 
metry  given  in  fig.  7  which  represents  a  particular  element  in 
the  band,  approximated  by  two  plane  triangles. 

Initially  a  coarse  scan  is  carried  out  just  to  detect,  but 
not  to  evaluate,  an  intersection  point.  This  is  done  checking 
for  each  sub-domain  whether  (ZX(I-l)  -  HT)/(ZX(I)  -  HT)  0. 

If  this  condition  is  true  a  root  is  detected  and  control  is 
tranferred  to  its  exact  evaluation.  The  value  of  the  root  is 
calculated  after  its  location,  either  in  the  first  (left)  or 
second  (right)  triangle  is  determined  (each  grid  rectangle  is  com¬ 
posed  of  two  triangles).  Equation  6  which  is  based  on  fig.  8 
(for  a  left  triangle)  or  eq.  7  which  is  based  on  fig.  9  (for  a 
right  triangle)  is  used  to  evaluate  the  penetration  point.  These 
equations  express  the  linear  interpolation  of  Cp  in  the  Fortran 
rotation  used  in  this  program. 

XS  =  X(I-l)  +  ABS( ZX( 1-1  )-HT)*(XM-X( 1-1  ) )/(ABS(HT-ZM)  + 


+  ABS ( Z  X ( I -1  ) -HT ) ) 


(6) 


Y 


7:  The  geometry  of  a  linearized  calibration 
surface  element  comprising  two  plane  tri¬ 
angles  and  its  intersection  with  a  plane 
normal  to  X,Y  along  Y  =  YG. 


abs(ht-  zm)> 


x(i-M 


xs- 


Fig.  8:  The  geometry  for  a  pe 
hand  (first)  triangle 
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Fig.  9:  The  geometry  for  a  penetration  through 
right  hand  (second)  triangle. 
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XS  =  XM  +  (X( I)-XM)*ABS(ZM-HT)/(ABS(ZM-HT)  + 

ABS(HT-Z(I)))  (7) 

Equations  6  and  7  are  invariant  to  the  slope  of  the  cali¬ 
bration  surface  i.e.,  the  slope  of  the  straight  lines  AB  and 
BC. 

Finally  it  should  be  pointed  out  that  when  two  penetration 
points  are  determined  the  values  of  their  abscissas,  X(I),  are 
returned.  When  only  a  single  penetration  point  is  detected  the 
abscissa  X(2)  will  be  returned  with  a  value  of  1000.0.  When 
no  penetration  point  is  determined  both  X ( I )  values  returned  will 
have  the  value  of  1000.0.  The  program  logic  is  designed  to 
recognize  these  messages. 

It  was  stated  earlier  that  PENPTS  is  used  to  determine  the 
closed  intersection  curves  projections  on  the  X,Y  plane.  In  fact 
not  the  curves  but  just  the  intersection  points  between  each  two 
of  them  are  required  (see  figure  3). 

To  compute  the  coordinates  of  these  points  the  subroutine 
"INTSCS"  is  used.  It  uses  PENPTS  as  a  subroutine. 

In  INTSCS  a  scanning  procedure  is  carried  out  from  a 
minimal  value  of  Y  (or  $)  YRIN  to  an  upper  value  of  YRUP,  or  a 


#>sss 


prescribed  10,000  times*  which  ever  comes  earlier.  The  sub¬ 
routine  scans  through  any  two  arbitrarily  chosen  closed  curve 
projections  to  find  their  intersection  points.  In  each  scan 
(I)  up  to  four  penetration  points  can  be  determined,  while  the 
penetration  points  of  the  previous  scan  (J)  are  memorized. 

Together,  eight  penetration  points  can  be  involved.  When  no 
intersection  point  exists  the  geometrical  situation  is  as  shown 
in  figure  11,  while  the  existance  of  an  intersection  point  is 
characterized  in  figure  12.  Subroutine  INTSCS  can  distinguish 
between  the  two  situations.  In  figures  11  and  12  the  case  of 
four  penetration  points  found  in  each  scan  are  shown.  The 
subroutine,  however  will  handle  any  possible  number  of  such  points, 
from  zero  to  four.  A  "no  penetration  points"  is  assigned  an 
abscissa  value  of  1000.0  by  PENPTS  as  explained  earlier. 

When  an  intersection  point  is  detected  its  evaluation  is 
based  on  the  geometry  shown  in  fig.  12.  The  eight  penetration 
points  have  the  following  coordinates: 


A0L( RAJL  ,  YRJ) 
BJR(RBJR.YRJ) 
B I L ( RB I L , YR ) 
AIR(RAIR.YR) 


BJL(RBJL.YRJ) 
AJR(RAJR.YRJ) 
A I L ( RAIL , YR) 

B I R ( RBIR.YR) 


*  The  10,000  scans  are  governed  by  a  program  constant  in  the 
line  "DO  140  ...."  and  can  be  arbitrarily  varied. 
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Fig.  11:  The  geometry  for  two  successive  scans 
when  no  intersection  point  exists. 


The  left  intersection  point  illustrated  in  fia.  12  is  the 
intersection  point  of  the  lines  aa  and  bb,  each  described  by 


its  equation: 

for  aa  Y  =  ( AA) X  +  BA  (8) 

for  bb  Y  =  ( AB ) X  +  BB  (9) 

The  constants  AA ,  BA,  AB  and  BB  are  given  in  eauations  10  to  13. 
AA  =  ( Y RJ- YR ) / ( RA JL  -  RAIL)  (10) 

BA  =  YR  -  AA*RAIL  (11) 

AB  =  (YRJ-YR)/( RBJL-RBIL)  (12) 

BB  =  YR  -  AB*RB I L  (13) 

The  coordinates  of  the  intersection  point  to  be  calculated  are 
X  =  (BB-BA)/(AA-AB)  (14) 

Y  =  X*AA  +  BA  (15) 


These  relations  are  true  for  a  left  hand  intersection  point. 
Analogous  equations  are  true  for  a  right  hand  intersection 
point.  In  this  algorithm  the  straight  lines  aa  and  bb  approxi¬ 
mate  the  curved  lines  connecting  AJL  with  AIL  and  BJL  with  BIL 
or  similar  lines  on  the  right  hand  side  of  figure  12.  The  error 
introduced  through  this  approximation  is  reduced  as  DAR  =  A* 
is  reduced. 

It  is  possible  that  an  intersection  point  is  identical  with 
AIL  and  BIL  or  AIR  and  BIR.  This  case  is  defined  as  "direct  hit 
The  program  is  designed  to  detect  such  a  direct  hit  and  evaluate 
the  corresponding  intersection  point  accordingly. 

The  above  algorithm  works  perfectly  as  long  as  the  two 
closed  Cp  =  const,  curves  are  far  from  being  tangent.  But  in 
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practice  a  situation  of  almost  tangent  curves  can  arise  when  X 
(or  a)  is  very  small.  In  this  case  the  preceeding  algorithm 
will  fail  and  must  be  replaced.  The  geometry  of  this  situa¬ 
tion  is  described  in  figure  13.  This  situation  is  identified 
by  INTSCS  and  the  intersection  points  are  then  evaluated  assuming 
that  they  are  intersections  of  two  circular  arcs.  When  the 
curves  Cp  =  const,  are  almost  circular  this  approximation  does 
not  lead  to  unacceptable  errors. 

It  was  stated  earlier  that  INTSCS  scans  from  a  minimum 
value  of  <j>  to  a  maximum  value  of  <f>  with  prescribed  steps 
a  $  •  as  shown  In  figs.  11  and  12.  This  direction  of  the  scan¬ 
ning  is  used  when  the  intersection  points  of  the  curves  Cpj 
and  Cpjj  of  fig.  5  are  evaluated. 

However,  in  course  of  the  reduction  of  the  measured  data, 
scans  in  the  direction  of  a  in  steps  of  Aa  are  also  neces¬ 
sary  to  evaluate  the  intersections  of  the  curves  CpI  and  Cpjjj 
or  Cpj  and  CpIV«  INTSCS  is  designed  to  carry  out  this  task 
as  well.  To  do  this  the  calling  statement  for  INTSCS  is  appro¬ 
priately  changed  as  will  be  explained  in  the  next  section.  To 
be  general  enough  INTSCS  is  not  written  in  terms  of  a  and  *  or 
X  and  Y  but  rather  in  terms  of  general  arguments.  The  best  way 
to  understand  INTSCS  is  to  compare  its  general  arguments  to  physi¬ 
cal  quantities  by  means  of  the  calling  statements.  In  fig.  14  the 
flow  diagram  in  INTSCS  is  given  in  terms  of  the  general  arguments. 

INTSCS  returns  to  the  main  program  (fig.  15)  the  (a,$) 
coordinates  of  two  intersection  points  between  the  Cp  =  const, 
curves  specified  in  the  calling  statement.  Let  us  now  follow 
the  way  in  which  the  main  program  is  designed  to  utilize  INTSCS 
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for  the  evaluation  of  o  and  $  of  the  velocity  vector,  as 
well  as  static  and  total  pressures. 
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Evaluation  of  the  Velocity  Vector  and  Pressures  from  the  Probe 
Si qnal s 

Fortran  program  VDR  was  written  to  evaluate  the  velocity 
vector  from  probe  measurements  of  pressures.  The  program  is 
shown  in  Fig.  15. 

At  line  1410  INTSCS  is  called  to  scan  curves  I  and  II 
for  possible  intersection  points.  Scanning  can  be  carried  out 
in  the  direction  of  the  ordinate  only,  with  the  calibration 
curve  matrices  compiled  exactly  as  shown  in  fig.  6.  This  limita¬ 
tion  is  imposed  by  the  way  PENPTS  is  constructed.  In  the  case 
of  the  intersection  points  between  I  and  II  ZA  and  ZB  are 
scanned  without  difficulties  in  the  direction  of  the  ordinate 
which  is  $  as  shown  in  fig.  12  and  returns  with  the  coor¬ 
dinates  of  the  two  first  intersection  points,  points  1  and  2 
of  Fig.  16.  They  are  ALF1,  PHI1  and  ALF2,  PHI2.  In  line 
1500  INTSCS  is  called  again  to  scan  curves  I  and  III.  Now 
scanning  has  to  be  carried  out  in  the  direction  of  the  abscissa, 
a  task  for  which  INTSCS  was  not  designed.  To  overcome  this 
problem  the  calibration  curve  matrices  are  used  in  a  transformed 
form  such  that  the  previous  abscissas  are  now  ordinates,  or¬ 
dinates  are  abscissas  and  the  internal  structures  of  the  Cppp 
arrays  are  accordingly  modified.  This  transformation  is  carried 
in  VDR  in  the  section  between  lines  320  and  450.  Comparison 
of  the  calling  line  1500  to  the  previous  calling  line 
1410  shows  very  clearly  how  the  various  arrays:  original  and  trans¬ 
formed,  are  used.  The  coordinates  of  points  3  and  4  of  fig.  16 


are  now  returned  to  the  main  program  of  VDR,  they  are  ALF3.PHI3 
and  ALF4.PHI4. 

We  are  seeking  a  single  intersection  point;  the  one  repre¬ 
senting  yaw  and  pitch  as  sensed  by  probe  A  in  positions  I  and  III 
and  probe  B  in  position  II.  But,  as  the  calibration  surfaces 
are  double  valued  we  are  now,  unfortunately,  in  the  possession  of 
four  points.  The  solution,  however,  is  physically  unique.  Only 
a  single  velocity  vector  exists  in  reality  and  its  yaw  and  pitch 
are  included  in  the  four  intersection  points  evaluated.  Were 
both  the  measurements  and  the  numerical  procedure  absolutely 
accurate,  two  of  the  four  points  would  have  been  identical.  But 
this  is  not  the  case  in  reality,  instead  of  a  sinole  point,  two 
points  close  to  each  other  will  be  detected.  Therefore  the 
average  coordinates  of  the  two  of  the  four  intersection  points 
returned  to  VDR  from  INTSCS  which  are  closest  to  each  other 
are  selected  as  the  measured  yaw  and  pitch  angles.  In  the  ex¬ 
ample  shown  in  fig.  16  this  will  be  the  point  dividing  the  dis¬ 
tance  between  points  2  and  3. 

The  calculation  is  now  at  the  point  at  which  a  and  q  are 
temporarily  known  (see  fig.  3).  Using  a  and  <j>  a  new 
CprpiV  (for  probe  A  in  position  IV)  is  computed  by  linear  interpola¬ 
tion  using  subroutine  INTPLT.  With  this  new  value  of  CprpIV  and 
with  PpIV  a  new  P$  is  computed.  If  this  new  value  of  P$  is  close 
enough  as  determined  by  EPSPSG  to  the  guessed  value,  or  to  the 
value  of  P$  in  the  previous  iteration,  the  data  reduction  for  the 
particular  point  in  question  is  terminated. 
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The  relative  difference  between  present  and  previous  P& 
is  compared  to  the  convergence  criterion  EPSPSG.  This  criterion 
is  evaluated  by  an  empirical  function  determined  to  give  best 
compromise  between  accuracy  of  results  and  computer  time  reauired 
until  convergence  is  achieved.  When  the  calculation  follows 
normal  routine  the  convergence  criterion  is  given  in  line 
1710  as  function  of  Pj  .  When  the  routine  for  almost  tangent 
curves  is  used  during  data  reduction  either  on  the  right  or  the 
left  side  a  different  empi rical  function  (lines  1720  or  1730) 
is  used. 

If  the  convergence  criterion  is  satisfied  results  are 
printed  out  and  data  for  a  new  measurement  point  is  read  for 
reduction.  If,  however  convergence  is  not  reached  a  new  static 
pressure  for  a  next  iteration  is  evaluated  (line  1780). 
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Convergence  and  accuracy 


The  convergence  of  the  present  iterative  procedure  is  not 
ensured  in  all  possible  cases  of  data  sets.  The  program,  however, 
is  adjusted  to  converge  in  most  of  the  cases.  Similarly  to  other 
iterative  computation  method  some  experience  is  required  to  achieve 
convergence  when  the  calculation  does  not  converge.  In  this  para¬ 
graph  the  principal  factors  affecting  convergence,  computation 
time  and  accuracy,  which  are  abviously  coupled,  are  pointed  out. 

The  convergence  and  accuracy  of  reduction  of  a  data  set:  Py, 
PPA,  PPB ,  PPC,  PPD,  PSIN  (or  analogous  set  in  the  experiment 
simulation  mode)  depends  on  the  following  factors,  which  can  be 
varied  by  the  user. 


1)  Coarseness  of  calibration  arrays  ZA,  ZB  and  their 
linearity.  With  coarser  arrays  convergence  problem 
will  increase  and  accuracy  of  results  decrease. 

2)  Probe  settings  YPB,  XPC,  XPD  (with  the  following  probe 
settings  being  fixed  and  not  variable  YPA  =  YPC  = 

YPD  =  0.0;  XPA  =  XPB  =  0.0).  The  optimum  setting  is  about 
+20°  to  25°.  Too  small  values  reduce  accuracy,  too  biq 
values  cause  convergence  problems  and  probe  tip  flow 
separati on . 

3)  RELXPS,  the  static  pressure  relaxation  factor.  The 
smaller  this  factor  the  safer  will  convergence  be 
achieved.  Computation  time,  however,  will  increase. 

4)  The  constant  5.0  in  line  1050  of  the  program.  The 
dimension  of  this  constant  is  kg/m2.  If  too  big 
results  can  be  lost  and  if  too  small  computer  time 
will  be  growing.  The  variation  of  this  constant  should 
be  coupled  with  an  appropriate  modification  of  the 
constant  in  line  1010  of  VDR  (item  8  in  this  list). 
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5)  The  constants  in  the  evaluation  of  EPSPSG  (lines  1710 
to  1730).  The  smaller  EPSPSG  the  more  accuarte  the 
results  in  expense  of  Increased  computation  time  and 
reduced  convergence  safety. 

6)  The  constant  scanning  step  DAR.  The  smaller  DAR  the 
more  accurate  the  results,  but  too  small  values  can 
cause  complete  loss  of  results.  Computation  time  in¬ 
creases  with  reduced  DAR. 

7)  The  constant  10  in  line  1360.  This  constant  governs 
the  number  of  Pg  corrections. 

8)  The  constant  1000  in  line  1010.  This  constant  governs 
the  number  of  scans. 


If  convergence  is  not  achieved  in  a  particular  case  varia¬ 
tion  of  each  or  of  a  combination  of  the  above  values  will  always 
enable  convergence. 


Evaluation  of  the  Computation  Time 

By  deleting  the  letters  CT  from  column  1  and  2  of  lines 
20  to  70,  1930  to  1960,  1980,  2000,  2010  the  actual  computation 
time  as  well  as  CPU  time  are  evaluated  and  printed  out.  The 
following  statement  prior  to  execution  is  required  in  this  case:* 

GLOBAL  T  SYSLIB  SSPLIB 

This  option  is  useful  for  adjusting  the  constants  affecting 
convergence  as  to  optimal  compromise  between  accuracy,  ease  of 
convergence  and  calculation  costs. 


*  When  the  program  is  run  on  Naval  Postgraduate  School  IBM  360 
sys  tern. 
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Input 


1)  Calibration  arrays  have  to  be  input  in  the  following  manner: 
first  values  of  NX  and  NNY  ,  second  values  of  X  in  rising 
order,  third  values  of  Y  in  rising  order  and  fourth  values 
of  Z  in  the  order  shown  in  Fig.  6.  Two  calibration  arrays 
are  read:  first  the  array  of  type  A  probe  and  second  the 
array  of  type  B  probe. 

2)  If  the  program  is  run  in  experiment  simulation  mode  PT 
and  PS  are  read  in  Kg/m2  and  subsequently  o  and  $  . 

3)  If  the  program  is  run  in  data  reduction  mode  the  measured 
values  of  PT,  PPA,  PPB,  PPC  and  PPD  are  read  and  then 
the  guess  of  static  pressure  PSIN.  All  are  read  in  Kg/m2. 
PSIN  has  to  be  lower  than  the  actually  existing  value  to 
ensure  convergence  of  the  calculation.  Too  low  a  value 
will  cause  waste  of  computer  time. 

Calibration  arrays  are  read  from  a  disk  space  on  which  they 
are  stored.  The  following  statement  prior  to  execution  is  re¬ 
quired  for  successful  reading: 

FILEDEF  02  DSK  NAME  XX 

here  NAME  is  the  name  of  the  file  on  which  the  calibration  data 
of  probes  A  and  B  is  stored  in  proper  order  and  format,  XX 
is  a  two  digit  number. 

Simulation  or  measurements  reduction  input  is  read  in  the 


normal  way  using  the  terminal  keyboard  or  punched  cards. 


a,  $,  Pg ,  Pj  and  Ma  are  printed  out  according  to  lines 
1820  to  1920.  This  output,  however  is  not  sufficient  when  the 
logic  of  the  computation  is  to  be  followed  either  to  examine  the 
execution  of  the  program  or  for  debugging.  Two  programs  with 
additional  output  are  given  in  appendices  2  and  3  and  can  be  used 
for  this  purpose.  The  first,  SWVDR  gives  short  additional 
output,  namely: 

1)  When  PSIN  is  successively  increased  automatically 
by  the  program  to  ensure  convergence,  the  values  of 
PSIN  are  printed  out. 

2)  1 1 T ,  FALF,  FPH I  and  PSN  are  printed  out  at  the 
end  of  each  iteration  prior  to  the  convergence  test. 

The  second  program  WVDR  prints  more  detailed  information.  All 

additional  WRITE  statements  in  this  program  are  numbered  with 

three  digit  numbers  starting  with  9. 


Conclusions 


The  following  conclusions  are  drawn  on  the  basis  of  experi¬ 
ence  gained  in  running  the  program  with  various  data  sets  and 
various  types  of  calibration  surfaces. 


1)  Optimal  probe  settings  for  ease  of  convergence  and 
high  accuracy  are:  YPB  =  25°,  XPC  =  25°,  XPD  =  -25°. 

It  is  therefore  suggested  that  a  B  type  probe  with 
25°  pitch  will  be  used,  and  that  the  A  type  probe 

be  rotated  to  +  25° . 

2)  Convergence  and  accuracy,  as  well  as  computer  time 
efficiency  are  improved  when  the  calibration  surfaces 
of  both  probes  are  not  flat  at  their  peaks  but  are 
rather  rounded.  It  is  therefore  suggested  that  a  new 
probe  tip  geometry  be  considered.  A  spherical  tip  with 
a  central  pressure  tap  is  recommended.  To  prevent 
damage  to  the  sensitive  transducer  located  behind  the 
pressure  tap,  and  in  order  to  improve  the  frequency 
response  of  the  probes  it  is  suggested  that  the  volume 
between  the  pressure  tap  face  and  the  transducer 

be  filled  with  an  appropriate  liquid  that  will  not 
affect  the  transducer  negatively.  In  this  case  the 
opening  of  the  pressure  tap  has  to  be  sealed  with  a 
very  thin  low  inertia  membrane. 

3)  It  is  probably  possible  to  modify  the  iterative  pro¬ 
cedure  such  that  safe  convergence  can  be  achieved  also 
when  using  the  scheme  of  Fig.  2.  If  this  can  be 
achieved  the  Kiel  probe  will  not  be  necessary.  An 
effort  in  this  direction  is  suggested. 
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APPENDIX  I 


LISTING  OF  VDR 
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07/08/74  id.A3.it 


FILE:  VO R 


FUR  T*4N  PI 


C  START  OF  VOR 
CT  JNTtCf  h  /)(6i 

CT  lAr>C*0 

CT  I A  tO*  L 

Ci^eRSlVl  aLf  lit)  i.’HI  .J  )>  X4  J  40), VM  AO)  , 
l/A»AC.'.Ci  ,A“  itOj  ,Yil40  I  ,2 i 140,40)  ,X  NX  140  I*  Y  AX  1 40  I  < 
27  \M  AC,  At  I  ,  A  >Xl  Av.  I  ,  t  iXl  40)  ,23X140,40) 

119  FC4mA’U1j1 

131  FOA/ATIF  7. 1 1 

130  FJ <M  AT | F  5. U 

120  FC  4  F AT  ( t-  d .  3  i 

1j2  Fir  MAT!  Fit.  3) 

135  f  GF'*  AT  (  F  13 .5  t 

13c  FJOVjTlfl/.jJ 

C  FEAO  PRO  ft  CALI -i-sM  ION  D1T\ 

Ac  A  J  (2,114)  .>A,NY 
DO  1<! 

121  REAClc, iJ)i  X  M 1 i 
C  J  122  1*1  , NY 

122  REM3(2tl3Cj  »Ml) 

CO  123  1*1,  i:< 

OJ  123  0*1 , 1Y 

123  R£AC(2'12J)  TM  I.Oi 
CO  221  l  *  1 » M 

221  REaOI  2t  130  x-li) 

CO  222  l*i,  IY 

222  FE  AC (2 , 1 )C  J  YHJI 

CO  223  1*1, NX 

CO  223  J  *  1 » \Y 

223  RE401  2, 120)  i->ll  ,  J> 

C  TR  AN  5F  CR  M  i  ALIF.-UTlon  MATRICES 
OC  224  1  =1  , N  V 

224  >1X111=  i4<  1 1 
00  225  1*1,  NX 

225  v*M  j>  *xa(  i) 

DO  226  1*  1, MY 
CC  226  J=i,M 

22  6  2AM  1  ,0)  *2A<  j  ,11 
00  227  1*1, NY 

227  XdX(I)=Ytli) 

JO  22 t  1*1, IX 

228  V  dx 1  I )*  X  E  1 1) 

OC  22S  1*1,  VY 
JO  2iV  J*  l, Mx 

229  28X11, J  1*2 -Mi, I) 

C  FRCe*.  StTYiNcS 

YPA-C.O 
VP  E*  55  . 

>p  a*  o.  o 

XPC*  20 . 

XPL*-2 J . 

C  PR3GRAF  Cu'ilA.MlS 
VRLn*-jJ  . 

»R1F*3 J. 

VR  IN*-  3  C  . 

X  FLC* - j J  • 

>RLP*30. 

>Rlfl*-jC. 

ICCFS*l 
NOCJPS*  1 
NO  S I M* 1 
11T«1 
iSCAF*l 
RCLXPS* J • 5 

C  EXPEM Ac  f  V  SIMULATION 
tO  1  FORVi'TlFli.tl 
193  Fl<C  15,  *01)  ->T 
FR11EI6,40l)  PT 
IF  I NCS  li.EO..,)  i  I  n  500 
ftt  AO  1 5 , 401 1  Pj 


IN  f  000  lo 
I  NTC0020 
I.1TC0C30 
I  1T0OJ40 
INTO  U1 50 
l  4TC006O 
INTO Jo  7 J 
IN  TCoCeJ 
INTCoCSO 
IMfOJIUO 
MTCCllO 
INTo^UJ 
INltltj) 
IN  TOO  140 
I N  TOol 50 
INTCCUC 
INTOJl  7J 
I N  TJJi Uu 
IMTlCISC 
IM0J2JJ 
1  NT  0U2 10 
INYC02  20 
INTO 02 30 
I N  TCC240 
INT002  30 
INTO 0260 
INTO 02  TO 
INT  002  <10 
INTCC2S0 
1NTCC2C0 
I. NT  00  3 10 
INTCC320 
INT  003  30 
I N  T C 0340 
1MTC0350 
INTUJ3t>0 
I.NTCC3  7C 
I NTOO  3dO 
NT  00390 
NTC04C0 
NT  004 10 
Nr0042J 
NTCC43J 
INTO  04  4  0 
I  1 T  0  o4  5  J 
IN  TCC460 
INTO 047j 
INTJjtRO 
INT 00 440 
I  NT 005 00 
INTCC510 
IN.C052J 
INTJ053J 
INT1C54C 
INT00550 
I  >N  T Oo5  o 0 
1NTC057J 
INT005OJ 
INTCC5SC 
NT 00600 
NTO 06 10 
M  TCCE20 
[N  i"  C  06  30 
INT00640 
, NTCC65C 
IN  T  006  60 
j NT00070 
NT  v  vC  Cv. 
IriT0069J 
INTCC7C0 
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FILE:  VdP 


FORTRAN  PI 


fcRITEtc.AOU  P'S 
*EAD(5,4JU  -L  *  5 
t>«lTC  If.  ,  VJl  |  ALFA 
AeA0(5,4Cll  Mil 

wall t (o, «ol J  “hi  1  . 

CALL  1  ,TPlT(aA,YA  ,1 A  ,  ALFA,  PH  II  ,NX,*IY,CPA) 

alk*vf\-ap: 

CALL  HT3LT  I*  a,Y  A.Zt.ALFC,  PHII.NA,  ‘JY,  CPC) 
AlFJ-AlFa-xPJ 

CALL  1  ir?LT  CaA.Y A,ZA,ALrO,PHII ,MX,  JY,CPJ) 
PH  la*  Fill  1-YPd 

CALL  I  llPLTUa.YC-  ,ic  .ALFA,  PHI  E, NX,  NY  ,CPB> 
FdYN=PT-P  i 
PPA=CF-«*PJYN»PS 
PP-3=CP  l*PJYw*-P3 
FFC*CF0«FjYN4P5 
PP0=CP JPPJYNFPS 
C  PEAO  Me  ASDPfc’l EM T 5  CAT  A 

IF  4MiS**.td.l  i  uC  TC  501 
‘OC  RtAdlt.li:)  I’M 
KKITHo.IjoJ  ?->A 
PE Ad (5,1 35 1  PP3 
WRITE!  c,  OfcJ  P°1 
FE20  (5,1  ill  r PC 
WRITE!  C,lJ„)  PFC 
RE*C(5,  12  5)  JP0 
WRITEU.lJu)  PP.J 
50 1  PE  Ad  (  5,  1 2  ;i  PS  l  I 
WRITE  lo, loo)  ?il'! 

IF (  ISCaN.o.. I )  TO  131 
18C  1F( ISCVUd  l.lddd)  JO  TO  305 

llT*l 
ICUPS*1 
ISCAN=1SC  VIM 
FSii\*PS  I  M5  .0 

181  PS*oilM 
F  SC*  F£1 n 

psf=ps 

C  CALCHLATtS  PCESSOAE  COEFFICIENTS 
30C  P0Y.4  =  PT-rS 

CPA*(PPA-Pi|/PC)YM 
CF0»( FF8-FSI/PJYN 
CPC*IPPC-PS) /P0»N 
CPC» (FPC-PS»/*CYN 

C  CORRECTS  Pi  ASSLMPTiCN  TC  EfiSJRE  CF  ARE  IN  CAI 
IFIFOr  lPi.c  J.i)  <j J  T3  301 
/APAX=-100UJ  . 

C  AlL  F>  *XJ-,y(ZA  ,2  AMAX  ,nx  ,NY  ) 

IF  1  CP  A  •  j  f  ...  .  'A  A  )  P3C  =  l.Jl*JS 
IF  IlFC. lT. ei^LA 1  BSC=1.J1*PS 
IFICPo.m  I.CAMax)  Psw*l.dl*P5 
2£  *tA-~  l  JjjO  . 

CALL  FAXJ-Y  w<2,lEvAX  ,'IX  ,  NY  I 

imcpp.jj.zuaxj  ps:=i.oi*ps 

2A**  IF*  l  J  J  JO  . 

CALL  F  l  NA  c  Y  ( t  i*  F,  f*  ,NX  ,  NY  1 
IFICF*  .eT.ZA'il'Il  PoC*L.NV'‘OS 
IFlCFL.Lr.A  J  MM  r-)C  =  U.',e<*PS 
lF(CP0.kT.2AAl  O  FiC*C.l<.«PS 
23Fr*UdJu. 

CALL  MiNAFV  (  le  U  t*>  »N  ,NX  ,  i\Y  I 
1FICP».L  I.C’.'III  P5C*C.V’J*P3 

if iPSc.E:.jsr  4  c  •  ~~  3 ji 

1C0PS=1  CCPSM 

IF  1  ICr,J  o  .  JT  •  Id  )  OJ  Tn  ISO 

FS*P5C 

p  S 1 3  P  ir. 

CO  Tn  jJO 

C  CALC  CL  A  ItS  T>J  AcF  6  F.li  ANCLES  FOF  C  'll' 


CALidRATIJl  RAflOE 


INT  JJ7  Id 
I  M0^/2d 
IN  TCC7J0 
INI Co74J 
INToO/5) 
INTOoFoJ 
INT  00/ 10 
l  N  TC0  7/1O 
INTCu  7S  J 
I  N  T  j  o  d  J 
I  (TCkeid 
1NTJJ1:.  J 
INldd  «3 J 
1NTC-P4 J 
INTO  J  J  *5  J 
1 . 1 T  C  L  e  u  d 
IN  CO  H7J 
I N  T  0  J  3  3  o 
IN  TC  C  S  SO 
IN  TJ.,4  0J 
INTdmO 
INTOOVid 
INTOJOJJ 
1NTCCS40 
INT  Go'iid 
INTOoSoJ 
1.NTCCS70 
INTJJNdd 
I  N  TJ09N J 
IN  I'C  ILCO 
INTO  10  10 
1NTJ1C? J 
IN  fold  id 
1  NT 0*040 
IN  TC  1C  50 

INTJI070 
INTUCHO 
INT01  0*5.) 
I  N  TO  1  lOvJ 
INTC1UJ 

IN  I'd l  led 

INTCllJO 
IN  fdll 40 
IN  Toll 50 
IN  TC 1 1  CO 
INTUll/O 
INTO l 130 
IN  Toil  U 
INT  JUOd 
IN  TC  1  e  Id 
INT  C It  iU 
I  NTu it JO 
INTCU4C 
I NTd l£  30 
IN  TOUjO 
INTO  IF  70 

iNtdlkOO 

I N  T  C 1 ?  SC 
INTdiiOO 
INTCiJU 
IN  TC 1 2 20 

INV01 £ iJ 

l  N  To  1 34  ) 
INT  J  13  id 
INVOUud 
IN  7C 1 3  id 
IN  V  L  1  3  6  J 
INTol 
INTO  1400 
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f lie 5  vub 


fchtuak  pi 


3Ci  call  nrs:si  yp;  ,ypi,y*:lc,y  f'jp.WI'i, 

1C  PA  ,CPc,  flLt-  ,  Jrli  ,  A  A,  r  4  U  \  , 

2X0.Yb,*£ ,A  AX.V-A , 2  AX ,Xtx, YK ,2fX,NX  ,NY , 1EPS  4 
ICPS1=(CPS 
/LH-PLM1  4 
4L  f  2  =  A  l  (-  I’l 
PhIl  =  Prl 1(14 
Ptl  i 2  *  Rn  I  u  I 

C  CAcClo.UCS  6  HI  AHuLtS  F3*  t.  •  III' 

CALL  1  1TSC.  ' '  ‘ 

1CPA 


2  >4  > 1 1 
IEPS2 


Irtw  —  fr  ii  r  It  h'lUkt  J  r  jn  I 

l  i  iTsc^ixPi.*0:,  >c.l\  .<puo,  r<n, 

tCPLifiil.ALl"  .X  (X  »V.!X  il  U  i 
i  Ki  a  ti4  t  tiA  t  lA  i  Vi  ilA  t^iY  ,  NX  t 
52= 1 t’5 


IEPS  4 


ALF  >**lF  t  l) 

ALF4=AcF( * ) 

Fhl3  =  F.U  (14 
Prtl4=FHH24 

SELECTS  TK  P-i  lPc-i  AlP  £  Phi  ViGLEi  ">'JT  DF  THE  F  JJK  VALLE  i 

FALEKlPlfK-Lr  Stiffs*  A  Lf-*  4—  AMJXl(ALFl»ELF2,ALF.a, 

1  *LF  <»  i-  A  A 1  t  A^F2  »-.lF  >  »A  Lr  X)  )  /  > 

EPK  =  (  (Pi-  11  +  'FU  *->h  l-tPhl  t  )-\  <4X11  PHI  It  PH  I  2,3  (1  A.  PHI  A  4 
1-AMI  ill  (  PM  1 1  *  P  H 1 2  ,  F.  1 1  J  «  3u  1  ♦  4  4  /  > 

CALCUL  AT  F  5  1E.1  Pi  K<‘.  v  J  »T  A  Oh  PVsiTlON  IV 
ALFCMFiLr-XPU 
CALL  IltlFLI(A.*,y-,Z 
PL  FCN=  f iLf-X  JC 


,ALFDi,FPrti,NX,NY,CPdJ 

PL  FCN=F  PLF-X  JC 

CALI  I  MUT  (  X  .4 ,  Y  P  ,  Z  A  ,ZLFCN»FPf  1,  MX  ,  NY,  CPC) 

psn=  (ppc  «•:  pl-pi’  u*c3c )/(  cpj-c  °c ) 

1  £  A  o  t  N  C  c  T  l  C  T  S 


CCN V£ A 1 1 A  - 
£PiPS= 


310 


303 

171 

17  2 
173 
1  74 


175 

CT 

CT 

CT 

CT 

450 
CT 

451 
CT 
CT 

3CS 
30  1 


EPj>P$  =  43it  l  Kn-P34  /PS) 

EP jPSC=  .COoCjooUoC;3<(PT**2.CI-.000  )003a*PT 
1KUFSI.il. 24  -P  5P  U=.JO  >  JO  30 JOJ 4* 1  FT  **2. j)-. 
IF  l  IFPS2  .E  ).l  4  i:  P  o3  SG  =  .  otJCu00L0O9*l  PT**  2.  C) 
IKEFSPS.oT  .tPSPou)  LC  TU  310 
GO  TC  3 C 2 
1 IT=  1  IT  ♦ 4 

IK  1  i  7.G  7.1  J)  jl  T;T  130 

Pi  =  PS«AELAPi*(  Pjli-P  j) 

FSC=PS 
PS  T  =  FS 
GO  TO  jJo 
V»K  I  T  t  (  j  ,  1  7  t  4  FALf 
FjXUTIU','  VtL'  CITY 
l  R  1 T  t  (o  ,1/  4  FPm 

FiJF  i*A  T  I  '  VuIClITY 

M  I  f'  c.  (  o  ,  1  7  .  4  3  o 

FCFKK*  i  T  AT  1C 
toft  1  TC  (  o  ,  1  7  l)  Pf 
FORMAT  (•  TJT.'.l  p 
4P*iLFT  1  (  l-T/P" 

WAHEU.l/tJ  i  X 

F3RM *T ( '  1  >Cf  IS'  «  F 15. :) 

CALL  IaLL’.nU! 

V  IF.T  1“=  (  A  15  4-  l  Ay  )/  /tdoo 
CPL=  iA(h4-  tif.v  4  /  4  63  10 
44IT£(u,45c4  VI  - 1  i  t 

FOF.'lT  (  •  V  i.^J-L  T  l"t  15  ' 

LAI  TE  (  *i  ,45  1)  CPU 
f  JR-fiTt  *  r  t  T \L  CPU  IS 
1  A6C=A( i 4 
I A  tC  =«  (  6  4 
CO  TP  19  3 
•  A  I  T C  (  j  ,  i  4  7) 

FC'V’ATl1  soACIeC  UJO  T  I  PP  S  '  4 
v»U  TC  19  J 
ENO 

SUt»0UT  HE  ("T  iCi  (  \P  ®  1,  p  15  2 , 4  "•  L  "*,  \PUP,A2IN 
lHTl,Hll  ,-l>1  ,Xl  |YI,11  , 

2  x  2  ,  V  fit.  ^  ,i  j,  y  3 ,  4  |4*4 ,  Y  *, ,  i.  •,  ,  C  i  ,  N  2 , 1  E  P  i ) 

Cl  MtNS  I  H  E  >  it  iJ  I ,  - t>  tl  i  j  4 ,  V'ujIIJIi 


JOUuOJc-’T 
vJOJOOJS^  °T 


YAV«  IS*  ,F  1  5.  2  4 

P1TC41  IS*  ,4-13.2  4 

RRESCUaC  LS»  ,  F15  .3  4 

EajC’E  Is'  ,r  13.  j  4 
.  J*=»0  ,2b>  71  4-1  .0  4*5  .0  4 


if  15.51 
,F  L  5.  54 


Xl(4o)  ,  Y  II  A  J  4  , 


INTO  IaU 
INTCi‘,20 
IN  TC  4*.  20 
IN  Td  It  *,J 
INTO  l  <*50 
IN  T  0  l  -‘*t  J 
lNTdli 7  J 
INTClAfcd 
ir.TUASd 
INTO  160 J 
1NTC1310 
INTO  1620 
I N  T  d  lr>  30 

14TCI  -J  *0 

INTO l 5 50 
|f|  TC15CC 
1NT01  5  7  J 

INT0156J 
IN TC  15  60 
IN  Vo  Ih  JJ 
INTCItlJ 
INTO  It  20 
1NT0U3J 
IN TC It  AO 
IN TClc5C 


INTOlttO 
I N  T  C  1 4  J  C 
I  NT'CU.dJ 
IN  VC  It  60 
I NT  Ll 7  CO 
INTC1710 
IN  TC1720 
INTO! 7j ) 
INTC17AJ 
INT C176J 
INTO  I  7 tO 
If)rCl7  7C 
1NV01  /  i>J 
INTO  1790 
INTCloOO 
INTO  It 1J 
IN  TCI t2f 
INTO!. 430 
INTCLE40 
INTCleEd 
INTO  I4'.o0 
INTClfc70 
INT  Jl  obO 
INT 01 6 90 


IN  TC  IS  CO 
1 N  V  0  1 6  U 
IN  •'014  20 


INT  Cits  ) 
I N  7  Cl  34  J 


IN  TCI  6  to 
INToiSoO 
IN’  CIS  70 
INI  013-0 
I  NT  Cl  V  SO 
1MTC2COO 
INT  C2Clo 


1 N  T  0Co20 
INT  O2o30 
1  N  T  02  d40 
INTC2050 
INT 02060 
I  N  T020  70 
IN  T C2  J  J'J 
INT  v2  0  6  5 
I  N TO t  ICO 
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FILE:  VJ3 


FIRTRAl 


1ZU  4C.4J)  .X2I40I  ,Y2l  40  J  ,  Z2  1*0  ♦<•*>  ,X3  <401  »Y3<40  »  . 
^  /j  I  <.0  >  -« C  )  #  *  *.  1  4J  it  (•ill  .CliZ,*I‘*0t4j) 

re  n  <  £  j  -  >  ioo. 

RtSl( i 

RE52  lll=9>>j. 

RE '2(21*51 73. 

IEPS*  1 
IS  1  =  1 
ISO*  1 
CAR=1. 

AA24*AAiK 
GO  I'?  152 

150  AR2R=AA2R *C4R 
Ar\2A  J*A12A 

152  AR*n  =  4-l2R-Atm 

Ai3P2*AA2A- JAP*. 

C  CHECK  FCP  lC*^.  CALI  cJaTllfl  RANuE 
JF  (  a*)**i.LT  .A-L  1 )  tl  TO  150 
IF|4PF2.LT.£rLC)  Li  TH  150 
K=  1 

C  CALCULATES  INITHL  P  E  1ST  ?.  AT  i  CN  POINTS 

CALL  rcAFTSUl  ,.\2  .<1  iYl  ,  Z1  ,HT1  ,ABU  .ARSES  1 
R AJL  =  A IRES  I  1) 

PA  »P*A  i«  S  tS  12  1 

CALL  PENPTS(;il,.s2  i*2  ,Y2.  li.  »HT2  ,AK'I2  ,  ARRES ) 

P3JL*  A  R r  tS  ( 1 ) 

PB jS=AkP  tS »2  ) 

C  CALCULATES  LLCCESIVt  “ENETRATl  CN  PCINTS 
CO  140  1=1.10000 
AR2P  *A i\2  x  ♦  JA y 
ARMi=iA2H- V3P1 
ASP2=AA2B-AKP2 

C  CHECK  FCR  LPPE*.  C*.Lt  EPhTI  CN  RANGE 
IF  I  ARM.GT.ARUP  l  Gi  TO  155 
IF  1  A  B  F2  .  u"  •  A  I'U  F  1  LJ  TO  193 

CALL  PE.MPTCINI.Au.aI.YI.ZI  ,nTl  »ARMl  .APKcS) 

PA  1 L=  A  S 1 1  S  (1) 

CALL  Ipc 'IP  i  01  iliNi  .*2  ,Y2,Z2  .HT2  »A R*2  ,  A  3PE  SI 
PBIL  =A JR£S ( 1 j 


TC  180 


RBI R  =  AA A  E  SI  2 ) 
iFIRAjL.fcJ.lOOJ. 1  I SL= 10 
IF  (PUL.60.1-JJ  ).  J  IS  L=  1 J 
IF1kBja.EC.  UOJ.  i  1  3fi  =  lo 
IF  l^JJH  .  kfl  i  Ijuj  .  I  l  j^=  10 
IFIPAIL.EC.I  >00. 1  l '  L  =  l  i 
IFIRFIL .E 1.1 JJU. )  1 SL  =  10 
IF  IPEIR.fc  J.iJOO  .  I  BA=10 

C  CHEC^Fns1  JHcCT  ’hit  *  Jl  ’lcFT  1  1TEA  SECTI  ON  »E  VALUA  TE  IF  RECUl  RED 
IFIISl.Kt.ii)  «C  TO  16a 
ISL*  1 
CC  TC  IdJ 

164  IF  1PA1L.Al.AA 1LI  LC  TC  180 
RESHK1=A  OIL 
PES2  <K>=AA2~’ 

IF  IK. fcO. 21  j  J  TJ  155 
K*  2 

GO  TC  13 j 

C  CHECK  FJR  l  .jTr  '*EJl  ATE  LEFT  INTERSECT  mu,  EVALUATE  IF  REO'JIRED 
180  IH  I.^L  .  It.  Lul  GO  Tj  JtO 
IS  l«l 

GO  TU  lii 

160  IF  ( l  =AJl-R*JI  )/ 1  :  i  IL-'RII  )  .ST.C.O)  GC  TO  ldJ 

AA*llr2rJ—  L  k  )  /  iPJJL~RAiLl 

04*4*  2  1- »3’=  <  lie 

A3*  (  A^^  F  J- J-'  J  >  \!  1  J  IL  I 

BJ  =  AA2*-**  A  >  i  L 

RESl  (K  1*1  ' ■*- ‘*.1  /I  Im-ARI 

FES2  lKJ*Aa*KtS14KI*EA 


luvcaio 

l NT  32 120 
INTl2U0 
1NT  02  1  40 
INW 150 
INTC21 tO 
INT021  n 

INTC21E0 
1M  T  021  .0 
INT  022  '0 
JNTC2210 
INT  J222  > 
I  NT 022  JO 
INT  J224J 
I  NT 02 2  30 
IN T  1.2*60 
INTC22  70 
IN  102*30 
|NTC*25C 
INT J£ J^O 
INTC2J  li 
INT C* i 20 
INTo23  Ji 
INTC2J‘.C 
INT 02 i  JO 
I  NTC23  'iO 
IN  T  1 2  3  70 
INT J2300 
1NT02JS0 
INT  G2**00 
1  NT  024  10 
1NTL2420 
1NTC2420 
1NTJ2440 
INTC245C 
1NTC2460 
IN  T 024 7  J 
IN  T  G  2  4  6  J 
I NTCZ490 
I.NTU25JC 
INT0*5  10 
1  NT  02  520 
INTC2520 
INT  C25*«0 
INT02R3J 
IN  T  C2  5  to 
INT  02  5  70 
iNT025y0 
INTC25S0 
INTC2600 
IN  T 126  1 C 
INT  02620 
INT026 >0 
INTC2140 
1.NT02650 
INT02660 
I  NT  0  2c 70 
INT 02630 
IN  TO* c 90 
INT027C0 
l NT  J2  7 10 
IN  TC2  720 
IN  T  J  *  7  20 
INT02  740 
INTC2  750 
INTJ2760 
INT0277J 
INTO.:  780 
I  NT  02  7  90 
IN  TC28C0 
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FILE !  VOS  FOOPUN  °1  N  P  S 

IFIK.CJ.2)  GO  TJ  135  INTC^glJ 

k*2  •  intuoimo 

C  CHECH  F  C  H  OlRcGT  hlT  Cl  RIGtiT  INT  ERS ECT I  IN, EVAL JAT E  IF  REO'IMEO  INTC2E30 

183  IF  (I  j>».  ME.  10)  uJ  TJ  lol  INT024.40 

IS  S=  1  I  N  T  02  3  'Jj 

uJTClM  1 N  T  02 1  tO 

U 1  IrIRiM.4E.Riri)  GO  TO  182  INT02070 

RES1  iR  !NT02tcG 

«Ey(M=iU2<  INTOcHJJ 

1FU.U.2)  GO  TO  15S.  I  NT  02  9  00 

M2  INTC2510 

GJ  TC  163  HVC292U 

C  CHECK  f'H  1  ITc  McOi  ME  RIGHT  I  NTER  SEC  TI  UN  »E  VALUi  TE  IF  RECUIREO  1NTC293J 

182  IF|I  iR.Kt.lO)  aT  TO  162  1NTC2S40 

I S  T3  1  INTJ2950 

GC  TC  1 6  J  IN  T 02960 

162  IF  (ItUOR-ReJ  U /(  Pi  M-RfiM  I  .GT  .0.0  1  GC  P  163  IN7C29J0 

AA*(  AR2RJ-A-.2R  )  /  4  U  jl-IA  M  I  INTJ2930 

fcA*AS2l-A  A''R  Ml  1)4102990 

AMUF2R  J-A12K)/  UUR-Roi  U  INTCOGOo 

PB*472<-AA*<- M  INTOjulO 

ft  SI  IF)  *lf  i-PAJ/  IM-Ae)  INTC3L20 

RES2lKU4.M«u.iHKU3A  1NTC3U30 

lFIK.EQ.2l  tP  TJ  1  Si  I.NT03040 

M2  JNTC34.SC 

C  RE-INITIATE  lNTOJJeO 

163  A«2fJ=AU1  IN  TOO 0 70 

RA JL*R  A  il  IN  TC30E0 

«AJR=MM  INT03J9U 

P6.iL=tdiL  IiNTOjiOO 

14C  R6j:<-R-m  INTO  j  1 10 

wRlTEto, 210)  1  NT  03 1  20 

210  FORMAT  l '  LOCP  1 4 J  THROUGH,  NQ  POINTS  F.1UNC*)  INTC3J20 

GJ  TC  15S  :  INTC314U 

195  If ( °fc? il  i).eJ.9G0C.C)  GO  TO  191  INTCJ150 

IF UESl 12) .cC. >000.0  )  Gw  rC  191  INTCJI6U 

IPUbf.21  ll.d.SCLC.C)  GJ  TO  19  1  INT0J170 

IF(«t<2<2).E'J.9GJj.Jl  G  I  TO  iVi  JNTOiliO 

GO  Tg  155  INTO JlSO 

C  Ey/ALUATIOI.  IE  I  -ITER  iter  IONS  WHEN  CLRVES  ARE  ALMOST  TANGEFT  ■  lNtOJCuO 

191  VGC=0.0  INTC3210 

I  EPS-=2  INTG3220 

CALL  FENPTi  I'12*N  1,aj,Yj,Z  J»NT1  » YGO»A  RRE  S )  IMTC32jJ 

PWiSI iPRtS (2 J-ARReS (1 ) )/2  IN TC 32 40 

CALL  Pt  J  2  *M  1  » Ai  , Y-*  ,2 4  ,nT?  «YGiJ»  4 R  AcS  1  INT0j25O 

R2*  US  <  A1  US  Ul- ATT  ESI  il  1/2  1NT0J260 

RL  J2ll)  =1  Ml**  »)-U2**2  MI  A1P2**2)  1/  I2.0*AT»2i  INT Oj270 

RE32I2  )  =  Rt:S2l  1)  INTJ32H0 

PtSl  IlJ  *SC«T  (US  I  Ml  **2  1  -I  7ES2  111**2  )) )  IN  TC2/-S0 

Rb  Sit  2>  *-.lu  31  t  l  J  1NT033C0 

155  RETURN  IN  T0 3  3 10 

EKJ  IN7CJ320 

SO  39  111  THE  P6NPT  SI  I X  »'•!  Y  *X  ,Y  *Z  ,HT  ,YG  ,XR)  INT03330 

0I*LK$1T4  Y  Ujj.ZXMOI.ZI  ,0,  40  1,  XI  4j)  .XRIlJ)  I.JT03J40 

11*1  INT033S0 

C  SEARCH  FOR  J  JF  LiME*.  y  LINE  INTO 3360 

GC  101  1*1, NY  INTCJ37C 

IF  (YMl.Gr.Yu)  ,C  TC  102  *  INT033SO 

101  CCNTINJE  l N  TO 3390 

102  j*l-l  IN  7 C 14 CO 

C  COaPjc  jG  a  i  FIR  ZE.-VJ  °1SS  INT0J410 

C  KtXT  4  L  I' '  c  i  ARE  ExECjrEO  II  FIRST  SEARCH  ONLY  INTOJ420 

YCM  YG-YI  J)  1  /  I Y  I  .♦!  l-Y  I  J)  1  INTC 34  30 

2x1  1  1*1  Zl  I  ,JM)-2  I  1, 0)  )••»:♦  2  I  1  ,J)  I  NTO  34  40 

lFIZXUl.tw.HTj  CC  TC  1 J 7  1 N 1  C J 4 SO 

13*2  INT034cO 

112  00  Kj  l=i'.K  INT  OoWo 

2XII  Ml  Zi  1  ,  JU  )-*  11,  J)  )*YC*ZI  i.J  1  INTC34SC 

Z'MIZII  ,  JUl-ZI  1-1,  J  1)  *(  Y..-YI  J)  1  /IYI  J*l)  -YI  J)  I  *Z(  l-l  ,4)  ’  IN  YOiViO 

IFIZxIIl.EUrM  wO  T)  luj  *  IN70J60J 
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file:  VJP 


FORTRAN  Pi 


N  I 
,  LTftilwO)  GC  TO  104 


if i i zxi  i-ii-hti/izxiii-mti 

10  A  CONTINUE 

C  T  C  lib 

C  CALCULATI"1  J  JF  oLNe  i'FATI  GNS  IN  EITHER  TRIANGLE 

104  . . .  . 


ri  ji  ) 


x*=x  i  i-i  i*ixi  n-xi  i-i i  i*t  yg-yi  J ))/{ vi  j*i»- 

lFiZF.tiJ.tiT )  UC  TC  105  .  , 

1FUZAI  I  1-HI1/1Z  i-HTl.tT.O.U)  GO  TC  1C5 
CO«AtS  (FT-i.^i  ♦.;  !>(  ZX  il-il-HT  i 
>S»XII-1I«-IXM-X1  i-1)  j*AbSI  ZXIi  -II  -HTl/OD 
CC  TC  110 

XS*X»<*  4X  11  l-AMjujdS  (Zv  -hT  1/  (  AfiSiZ'-FT  KA9SIZX1  1 1-FT  J) 
G J  TO  IK 
C  LOGICAL  Ft  MtTRATiJ'iS  ACCUMULATION 


105 


1C  7 

>*m*xm 

i*i 

GG  TP  113 

10E 

XKiiPKXl  (j 

go  n  111 

110 

XR  1  1R  )*XS 

111 

IF11R.E3.2)  uO 

TO 

114 

113 

1  R*Z 
ie-141 

GO  in  IK 

116 

IFilP.LC.l)  GC 

TC 

ii7 

XR 121*1000.0 

GC1  TC  1 1 4 

117 

XR  1 1 1 *1000  .0 
XR12>=ICCC.C 

114 

RETURN 

END 

SU^ROLT  i  Jt  I'Lr?L1(X,YfZ.XG,yGf-\X,NY.ZRES> 

C  LINEAR  iNTcRFCLATnN  CN  CaLiBR AT  1  .'IN  SURFACE  TO  EVALUATE 
C  CP  VALLE  aT  Ci.-i.rtJl  NaTcS  XG.Yo.THt  RESULT  RETURNeJ  iS  ZREi 
oi Pension  a (4J) , y( 4oi»zi 4u,4C) 

C  SEARCm  fCR  I  CE  LeFT  X  LINE 
DO  1  IC*1,MX 
IFlXilCl.CE.Aul  Cl  TCI  2 
1  CCNTIME 

Z  l*  1C-  l 

C  SfcAPCf  FCR  J  CF  LCWEF  V  LINE 
OJ  3  JC*  1  »  <Y 
If IYIJC l.Ge.YG)  3)  TO  4 

3  CONTINUE 

4  J“X- 1 

«L*lYU*ll-YlJ))/IX< 1*1 )-X(I)i 
BL*Y  IJ1  -Al<-'X  1 1 1 
VCR«<AL*XG)  ML 
IF  iYu.CC. YCK J  31  T1  7 
IF  I  Yu  «  j  1 .  YC  rt  j  Vj  u  lu  3 
C  RESULTING  Z  IS  I  4  LJRRER  TRIANGLE 

xj-xuhi 
V3*  Y 1 J) 

73* Z  iI*l,J  1 
GO  TC  0 

C  RASUL  T INC  Z  15  I  J  UPPER  TRIANGLE 

5  X3*X ill 
>3*V< J»l 1 

C  CENEFAl  CALCUlATICNslCCC  FCR  EOT  F  TRIANGLES 
«  A 1  *  X  i  I  1 

Y 1=Y  i  J  I 
Zl-ZIl »J> 

X2  =  XlIM) 
t2*y  ij*i  i 

ip*{ 11  1—  2  J  »  * < Yl-Y2)-iZ l-Z2)*<  YI-Y 3)1/1 <X1-a31* 1YI-Y2 J-1XI-X2 1 
I* ( Y1 -Y3  1  I 

BP*I  I ll  —  LZ 1  - «P* I X1-X21 1Z1YI-Y2 1 
CP  =  ZI-AP*aI-<P*»'  1 
7RtS*AP*Xu»ev*y G«CP 
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INT  0  35  10 
INT035Z0 
IN  TCJ53C 
INT  0  3  541) 
1NT0J950 
IN  TC  3st 0 
INTJ157J 
1NT0356C 
IN  T  J 3540 
IN  TO Jo JO 
IN  TCJt  10 
INTO  Jo 20 
INTOjt JO 
IN  T CJt4C 
INT0JO50 
INTCJf.60 
IN  TO  3e  7u 
INTO jo  30 
INT03690 
INT  JJ  / 30 
I  NTOj  710 
INT  C J  720 
INTO J7 30 
INTC3740 
INI  03750 
I  NT  0 J l oO 
HTC37  70 
IN  T0J7  30 
I  N  TO 37  NO 
INTO 3800 
iNTOjolo 
IMTC3S20 
Wf  0  3e30 
INTO J84J 
IN  TC3E5C 
INT  03  d60 
INT03870 
INTCiueO 

INTJJobO 
INTC35C3 
INTO  J9  10 
I  NT Cj5Z0 
I N  TC  j  5  30 
INTC3S40 
INTO J9 50 
IN1C3VO0 
INT 039  70 
IN  TO  3930 
INYC3SS0 
I N i Oho  00 
IN  TC4CIC 
IN  TC4020 
INT040JO 
INTC4C40 
I  NT  040 50 
INT  04060 
INT  040  70 
I  NT 04080 
INTC4CS0 
INTC41C0 
INT  041 10 
INTC412C 
IN  TO  41 30 
INT0H14J 
INTC4150 
I N  I  04 1 b  J 
INT0417) 
1NT04  1  eJ 
I  NT  04 1 -*0 
1NTC42C0 


1 


J 


Flit:  VUP 


fortran  ®i 


GJ  TC  IOC 

C  PESULTMu  Z  IS  )  4  Jivur.j  LMC 

f  s  »ui  i  n  ,j*i j-zi i,  j  j)*ix«-x  i  m/(xi  ini-xt  tn*zi  i  *  j  » 

10 C  ec foRN 
ENC 

fteULUNt  X  4  X  A  R  V  I Z  »  2  *<  AX  »  NX  •  NY  4 
C  EVALUATION  OF  >\xHlH  CP  Jf  CALIBRATION  SURFACE  NaTRIX 
CIFEFSlCA  ziAJ.AOj 
00  1  1*1,  IX 
03  1  J-»l,  4Y 

1  IF  Utl  ,  J)  .  ot.Z-'AA  4  Z*AX»Zl  I,  J  I 

RETURN 

lu-OlUNE  y  1  Nm 3  Y  IZ  »ZHIN,NX,NY  ) 

C  EVALU4  r  I  "IN  OF  MUI-IJI*  OP  CF  C-.  Ll  tik  ATI  ON  SURFACE  UTR1X 

I  li  L  A  C  I  '  fc.  I  4  .  I 


C  1XENSI0N  a14u,ij4 
0'J  1  I  -  1  ,  A  X 
03  1  J  =  It'l> 

IF  (Z  1 1 1  J  4  •  Lw.  ZM IN  4 

PE  TUAN 

ENC 

ENU  CF  VCR 


ZMi  4=ZI  I,  J  1 


INTC42U 
INTJ4220 
IN  TC42  30 
INT0424J 
iNTCtZSO 
INTCAitU 
INTJ42/J 
l.NTCAZtC 
IN  I  042  *0 
iNTCtJUO 
IN  TC431J 
IN  TJt  320 
I  NT  J  •*  i  >0 
I  NT  0<34O 
IN'.' JAi  jJ 
IN7C4KJ 
iNi C A  3 Jj 
In T  ovo  ;j 
I II T  C  A3  SO 
IN  7  0  *400 
IN  7 OVA  U 

In  t  c aazo 

INTjVA JJ 


ai»4  UIUIUIIII  » JIXXXXXAAXXXXXXXXXAXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXI 


APPENDIX  2 


LISTING  OF  SWVDR 


CJ/Ct/JS  It. 45. 15 

FI  It !  SrfVCR  FCRT  SAN  Pi  N  P 

C  STAR  1  3F  S« VJrt 

Cl^ENS!  14  AUF(  lo),  3  Hl(  10),  XM  40  ,  VA<  4C)  . 

1ZA140  ,40  J  ,<d  (4  0  J  ,  Y214J),  L  H43,-«J  ),XAX(4J  It  YAX(  40  4  ,• 
2ZAM  4c»  4C)  ,  X*  A(  -.O  ,  YMXUUI  ,2  3X140  ,40  4 

119  f CF“ AT  1215  I 

131  FuaVMC  F7.il 
13C  f.-!K4AT(f  5.  i) 

120  FCi'*A_<  F6.3  J 

132  FCH«AT(F  1,).  54 

135  FCF‘-1  ATI  F  li  . 3  J 

136  FCt..<*AT(Fl  7.54 

£  READ  PPJHE  CU.198ATI  /l  DATA 
REAC(2,ilii  IX,  4Y 
DO  121  1*1.  IX 
RFAQ(2,l30)  XU  1  ) 

121  Ci-NTU-JE 

01  122  1*1, Nr 
RtADU,  120  )  YA( I  ) 

122  CCNTlNJt 

00  122  1 *1 ,3A 

CO  123  J*  1,  4  Y 
FE  AO (2 , 1 20  I  24(1, oi 
122  COlOlMUc 

CO  221  1*1, MX 
REAJ(2, 120)  Add) 

221  CCMIMJF 

CO  222  1*1 , NY 
READ ( 2, 12CI  Y91I1 

222  CJNTirtiE 

OC  223  1=1, NX 
00  222  J*  l ,  4  Y 
FE Al (2 , 120/  2 E ( 1 , J 1 

223  CONTINUE 

C  TRANSFORM  C  AL I bR  N  T1  Jii  MATRICES 
CO  224  1=1, NY 

224  >1M()*UUI 
00  x25  1*1, NX 

225  YAXm=XA(iJ 
00  226  1*1, NY 
on  2  2b  0*1,  4  x 

226  2AX11  ,  Jl  *Z*UJ,1 1 
OJ  227  I  *1  ,!W 

227  xax(l)=Yfc(l) 

CO  22b  i*l , N X 

22  E  Yd>m*XPIl) 

CO  225  1*1, 4 Y 
00  229  J*1 ,  rX 
225  Z<)X(  J,JJ*Z?(  J,i) 

C  FRCBE  SE.  rfi-ics 
YPA*C. 0 
YP8*  55 • 

XPA=0. J 
XPC*  2C. 

XPC* *20. 

C  FRCOFAF  CONSTANTS 
Y»L0*-3C. 

YRU°*3 j . 

YRIN=-30. 

XRLC1 -*-30. 

XRUP*3J. 

>R  I N  *-  3  0  . 

ICCPS*  1 
NCCt  FS*2 
>10  SIM*  2 
I1T*  1 
1 SCAN*1 
REl.  >PE=  C  .  3 

C  EXPERIMENT  il-ULST  |U 

153  CTNUN'Ic 

401  FJFM AT I F 15 .4 ) 


NTCuCIO 
I N  TO  J JZO 
1NTCCC3C 
INTO  JO 40 
-  N  T 00050 
iNTCJCbJ 
I  NT 00 J  70 
iNTCCCfcO 
INTO J 040 
1  f4T'J  j  1 00 
IN  TOO  1 10 

In  T  jo  1 20 

1 N  T  Jo  1  JO 
INTO J 140 
I  iJTJJl  30 
INTCOltO 
IN  i  OC  l  70 
1 14  TOOl  d  ) 
IN TC Cl  50 
IN  f  002  00 
INU021J 
IN  TO 0  2 20 
1.IT0J2  iJ 
1NT0J24J 
IN  TOO  2  50 
INTJ026J 
1NTU02  40 
INT  00260 
INT0v>25  J 
INTCC3CG 
INTO  Jo  10 
IN T  00320 
INTCJJaO 
IMT0o34o 
1NTCC35C 
NT  00 jtJ 
NT  003  70 
NTC02UJ 
NTCC35J 
N  TOO-,  JO 
1NTCC41C 
IN  To042  J 
I N  TOO'*  30 
I4TCC440 
INT0045J 
NT  00460 
INT 004  70 
i NT  Oo4bO 
I  NTCi/4  90 
INTC05C0 
INTjo510 
INT  C  0520 
INTO  03  2-4 
1  N  TO j5  4  J 
1NTCC550 
INTJoNoO 
INT005  IJ 
INT0-5eo 
INT  0 0550 
IN  T  COt CC 
1NTCC610 
INT0J620 
INTCCt3C 
INTO  J640 
INT0O650 
INTCbctO 
INTOJoTO 
INT0J61J 

In toovvj 

I  NT 007 00 
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FILE  i  SMVON 


FCRTRa.N  PI 


•!!!*#  V 


i,2A,ALFA  ,PHll  ,NX,NY,CP4) 


WRITE ( a, 4J 1)  PT 
lFIKSiP.tC.2)  JC  TO  500 
RE  AO ( 5,  4l  l)  PS 
W  R ITE  16*  -»0 1 1  PS 
REA0!5,4Cli  ALi-A 
WR  ITEI  c,  4v.  1)  alF»* 

REAC  !5,40l  )  PHI 
kRlltlc.^l)  Prtii 
CALL  INTFl T (*t,r 
ALf C*ALF A-aPC 

CALL  li\ULTCXrt,»i.,i:fi  ,ALFC#PHII  , NX, NY, CPC) 
FL  FC*  AL  FA-APC 

CALL  IMFL7  txA.Y  /  ,2  A  ,aLF C,  FH  1,‘IX,  NY,  CPC) 
PF.18*PH1I-YP3 

Call  intclt  uh.y  f,2B,alfa,phib,nx,ny,c»0) 
P0YN*PT-F5 
PPA*CPA’*PJYN»PS 
FFt*CFj*fjYN+PS 
PPC«CPC*PJY  >i*-PS 
FFC*CP0  4PCY'I*-US 
READ  MEtSUPcPtNT J  CATA 

IF  INOSH.t  0.  i)  S  I  TO  501 
REAU5,iJ5)  P° A 
8RI  lEi't  ,136)  PPA 
READ) 5,133)  PPA 
t>  A  IT  £  !<j  ,  1  jo  )  PPO 
READ (5,125)  PPC 
WRITE!  0,  Uoi  PPC 
RE  AO (5 , 1 35  )  PPO 
WRITE! c,  lie)  PP u 
Ft ACI5, 1J3)  3  S  IN 
RRi T 1 16 ,  lib)  PSiN 
IF! lsrAN.EQ.tJ  JO  TO  181 
IFliSCAN.CT . iOOJ  i  oO  TO  305 
I!  1=1 
ICCPS-1 
I  SC AN  =  1 S  C AN*1 
PSlN*PalN*5.C 
W  A  IT  £  (  i  ,  )2c)  PS1.J 
CCMULE 
PS-PSlN 
FSC-PS1N 

C  CALCULATtf  PRESSURE 
300  CONTINUE 

PO  YN*P I -PS 
CPA* !PPA-PS)/pGYf 
CFE«(FPE-F£) /PCYN 
CPC*  i  PPC-P  S  )  /°UY’i 
CPt*  !  PP  C-:,S  )  /  a  CY ,'i 
C  CORRECT*  P 5  ASSLpFTIlN 
1FIMJC0PS. £0.2)  JO 
2ANnX*-l JOoJ . 

CALL  MAX4rlY!2A,2AHX,NX,NY) 

If  1CPA.CT .2*4  AX )  =SC*i.Ji*PS 
IF  1C  FC. u 1.  t  *"i  AX  )  F SC  *  l  .0  1  *  ri> 
lFlCPO.oT.ZA-lAx)  PSC«1.01*PS 
2BHX*-1J  )0d  . 

CALL  NAXApY(2b,2eN4X  ,NX,NY) 
IHCPR.CT.CS-iaa)  »SC*1.U1*PS 

2ANI N*iUOoO. 

CALL  Ml  N  A  i»Y  124, /API  i,  , NX,  NY) 

1MCPA.L  T  ,2A  1  IN  )  OSC  =  O.VN*PS 
IF  jCFC.LT.2JHN)  FSC*U.<»9*PS 
I  F l CPP .L  1 .2  4  4  I  I )  P  3C  *  C.99*  PS 
2B«lN*i )000 . 

CALL  MI.4APYI21  ,2dMN  ,M  ,NY  ) 
iFICPP.LT.23'1IN)  PSC«0.<)V*PS 
IF  1PSC.EW.PST )  LC  TC  30i 


500 


501 


180 


1C1 


COEFFICIENTS 


TC 

TO 


ENSURE 
JO  1 


CF  ARE  IN  CFL  i SR AT  ION  RANGE 


1NTCC710 
lNTOJ/20 
INTCCJ3C 
IN  T  00740 
I  N  TO J 750 
IN TC 0760 
1NT00770 
1NTCC7B0 
IN  !' JO 
IntcosO) 

IN  TC0810 
1UTCC82J 
I  NTJuT  tj 
1  M  TCCfc  '.C 
INi'd'JdiO 
INI  CuobO 
IN  I  Cot  70 
NTOJddJ 
NTCC8S0 
IN  T00900 
INT0J9I0 
INTCCS20 
IN  Tud9  30 
INT00940 
INTojySO 
I  NT  00960 
1NTCCS70 
NTCC580 
INT0J990 
IMTC1CC0 
IN  TO  10 10 
INTJ1 020 
1NTCIC30 
N  TO  1040 
INTC1C50 
INTGlObO 
I\Tvil070 
INTCICeo 
INTC109J 
NTOilOJ 
IN  TCI 110 
INT01120 
IN  TO  IliO 
INTCll 40 
INTO  1 150 
NTCllfcC 
INT011  70 
IN  701 1 30 
INTC11S0 

Into 12 oo 

INT012 10 
INTC1220 
NT  012  30 
IN  TC  12  40 
IN  VO  12  50 
INTC120J 
I N  TO  12  70 
INT  0 1280 
INTC12S0 
MTCUCO 
NT  0 10 10 
INTC132C 
IN  TO  1 3  30 
IN  T(J1  j40 
NTC13S0 
I NTU 1 J60 
NT01370 
NTOlOdd 
NT  01 j50 
INTC14C0 
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me»  swvcr 


FflRTRAN  »l 


o7  TO  lt)0 

pfr*p§c 

bu  TC  j JO 


101  (JonTIKCc 
C  CALCULATE*  it. 'I 
CALL  lwT>C3 


AcF  6  PHI  VlGl.cS  PC'!  '  I  '  t.  'll* 

,  __  _  JlrPA.Y'M.YRLT.YRJP.YwlN, 
lCPA  1CP3  f  ,  PH  I  ,  X  A  |V  A  f  i  a  , 

2XB.YR.Z A ,*A  X, YAA.ZAx.AdX.YHX .ZeX.NX.NY) 

ttfakthi 

PHli«PFi<  11 

c  calcljlHfs'1  l  i^ALf  C  p.ll  ANaLt  S  FOR  '  1  •  F.  'Ill' 

CALL  InTSCS  UPA.XPC.XRLO.XRIJP.XRIN, 

1CP  A  tCPw  (Pi  1  1  4  cF  .YbX.YAX.ZAX, 

ZXAX.YAx.ZTX.XA,  I  A  ,ZA  t  A*.  ,  YA  .ZA.NY.  AXi 
ALFa«ALC  111 
ALF  4«*»lF  121 
FHl3*Pnt 111 

pm  Whim 

C  SELECT*  THE  PTIPfc*  AtF  &  PHI  ANGLES  CUT  CF  THE  hOIR  VALLES 
f  ALF*  I  l  AlFl  *ilF<.  YALFjYiLFH  l-AMAXK  ALF  li  ALF2,  AlF3, 
1AV-4)*’A*|.*1  ( iLFl  ,AlF2  ,iLF3  .ALF411/2 
FPPl-l!  PHI  1+Pril2+P  II  3+PHi  4) -A^  Axil  PH  1 1 ,  PHI  2  ./Mi  J.PM41 
1  -API  M  (  PH  4  tPHlcF.1 4  A,  Phi*  1 1/2 
C  CALCULATES  JE«  »*  Fk.M  0'*T\  CF  PUjIT 1  OH  D  AMD  COXRfcCTS  PT 
ALFCN«r ALF-X’G 

CALL  1  ■nFsTIX‘*,Y*.,ZA1ACFDHlFfHI,NX,NY,CPCl 
AL  FCH«FAth-XPC 

CALL  1  i T o|_T  IX  A ,  /  A  ,Z  A  ,  A  LFC  .*l-»  P  PH  l,  IX.MY.CPCI 
PSN-IPPC+C  PJ-PPJ*CFC1  /(CFu-CPCl 
907  WRITE! a. 4  1  111  IT, FAtF, FPHl . P*.N 
91 1  FURPATf  ITfewaTlCN  VALUES1,  i»,  10X , 2F6  .2, F 10 . 21 
C  CflNVEKGENCc  f c jTS 

EPSPS*  AE  S I ( p  3i-P j 1/P  S 1 

-  -  -  ■ -•JT**2-a»-.o?>i)tiiiJoa*PT 

0  J'J0J8*PT 

C-OOCCE*PT 


310 


m 

172 

173 
114 


175 

3CS 
30  7 


GO  TO  jJJ 
CONTINUE 

111*11  1+1 

IMIlT.tT.iJ#  GT  T 1  180 
F  S* F  S+  * tiLX PS*  IPSN-PS  J 
PSC*PS 
FST-PS 
GO  TC  30C 
WRITE!*.  I/ll  FALF 

PC  At  I  •  0  •  »  '  VtLLCiTr  YAW  IS  •  ,F13.2) 

KUTclu.  1/21  F  Pm  1 

FCFMAT1*  VtLJClTY  •■»  ITCH  IS'  .F15.21 
NRiTt(6«l/Jl  PSN 

FORMAT!  1  ;)AilC  PRcSSImc  IS*  .F15.3J 
WRITER. L7.1  FT 

FORMAT!*  7  J  TbL  PnESSLRE  IS'  .Fla.il 
AM-SCFT  (  (  <.’T/P 3,' 1**0.21)57  1 4-1. Cl  *5.0 1 
hftITt  to .1  lot  AM 
FORMAT (  *  •I  TCH  I40F.9ER 
CO  T*»  14* 

WRi  Tt (o ,307  1 

-FORMAT!  1  jCAMMEU  1000  TIMES'I 
oU  TC  LOO 
ENU 

SU8RCUTHE  I  ITSLa  (  4RP  1,  4RP  2,  AR  LC,  1RUP*A3IN» 

1HT  l  »hT2  ■  b  tS  1  ,  rtjti.U.Y  1.  21  » 
2X2,Y*,2<,A2,Y3,22,X4,Y4,Z/t  .Ml  ,N2) 

Cl RfcNS  ICY  wC3l(lJi,R  ES2!  10  I, A’ RES!  101.XU401  ,Yl(  401  , 
1Z1  (40,40  ,X2  (401  ,Y2  (401  ,Z2  (40,tJ),X  3  <401,  Y  3!  40  l, 


IS'  »F l 5. 5 1 


MTC1410 
NTO  l*»  20 
1MTC143C 
'1  TO  14  40 
NT0143J 
hUlKO 
1NTJ147J 

In  to  1-.  50 

Nib  14*0 
NT01500 
IN  TC  15  iu 
1MT01520 
NTO1530 
I NTC154C 
INTO  I \t 50 
INTO  15  00 
TIT C  15  70 
NlJladJ 
IMTC15S0 
IN VOlc wO 
I  NT  016  10 
INK  1620 
NTJlbaa 
NT3lc40 
|NTCl65J 
I  NT JittO 
1NTLUT0 

IN  T  j 168J 

In  KIOTO 
INTO  1/00 
INTOi/lJ 
,  N  1 J  17  20 
IN1C1720 
NTOl/40 
IN  TCI  7 50 
INTO  1 7  wJ 
IN  TCI  7 70 
1MTC1720 
INTJlf-iO 
NT  C 1 1  CO 
MT31S  U 
NTOldZO 
M  TO  Id  30 
NTO  10  40 
NTvHoaJ 
IN  TOlHtJ 
I  NT  J 1 6  /0 
lNT0l«8o 
IMTCleSO 
INTO 19 00 
IMTCIS10 
1N  +  01R2O 
I  NTOiSiJ 
IMTCIS4J 
NTJ195J 
INTClStC 
IN  TO  14  /O 

iMToisao 
N701SS0 
INT020 JO 
NT 020  10 
,  NT H2U 20 
INTO 20  JO 
I N  TC 20  40 
NT02C50 
lNT020o0 
IN Tt2C  7C 
IN7020FO 
NT 02 040 
NTC21C0 
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Flits  SPVCP 


F3RTR4N  PI 


2.1  M  4C,iO)  ,X4(40J  (74(401  ,24(40*60) 

AES  1  ( 1 1  *  »  jj  w/ . 

Be  SI  (?)*4t)oi). 

Rti2(  II  =SCOC. 

RES2 (2  J«9.»00 . 

1SL«1 

isr*  i 

UA  P  *  1 . 

AKiR*A-<I  N 
bC  TC  152 

ISC  Ah?R*AA2«»C4K 
4R2RJ*  \k23 

152  AK'n*AR2R--RPl 
ARN2*a,«2P-4RP2 

C  CHECK  FOR  L'J«e.<  CALIBRATION  RANOfc 
|MJSn.LT.«KLOI  tC  TC  150 
IF  (  APF2.  L  !•  AKLCI  oC  1C  150 

C  CALCULATES  INlTUL  F  LLcT  F  AT  1 CN  PCIFTS 

CALL  PCNPTSl  a, .<2  ,Xl  «Y1  ,i  l,hTl,AR«l,ARR£Sl 
PAJL* APR  ES (11 

RA,iR=AAK£Sl2)  _  . 

CALL  PENPTS  ('a,N*,X2»Y2,2  2,HT2,AR"t2  tARRES) 

KB  «.l=A<"  PtS  4 1  I 
RB jP=AkRES(2) 

C  CALCULATES  SuCCESIVE  PENETRATION  POINTS 
LG  LAC  1*1,10000 
AP*R=AA2P«"J4* 

ARM1»A«2S-APP1 
ARP2  =  AH2H-r.rvP2 

C  CHECK  I  OR  UPPER  CALIBRATION  RANGE 
IFiAFf  i.LT.^UP  J  CT  TC  ISO 
|F(ARE2.jT.i>KUE)  O:  TC  19b  „  . , 

CALL  PL  'IPT  j  (  a,  N2t  XI,  Kit  i  l ,H T l  ,ARMl , ARAE S) 

RA1 t*Ar  Fc  j (1 1 

CALL*PtN?Tsf,H,N2iX2.Y2,2  2 ,HT2 ,ARM2 ,ARRE  S) 

PHILIPPES  III 
RB I R=ARR t  S  (  2 1 
IF  (PUL  .E  J. 1000  .  4  IS  L1 10 
IF (RtJL.EC.l  JJJ. J  IS  1-10 
IF  (RB  JR  .fcO.  1000.  I  (S3=lJ 
IF  (MJR.EC.  UJJ  .  i  ICS=io 
IFCP.AIL.EC.  10CU.  I  1  51*10 
1FIREIL .tO. luoj.  1  I SL=  10 
IF  ( F 1 1  R.  t  C  .  10  )0 . 1  13w* 10 
IF (R  All . 6  5. luOC. 1  IaB  =  10 
C  CHECK  FCF  OIHECT  MT  H  eEFt 
iFIlSL.Nt.lJJ  ^  TC  164 
ISL*  1 
CC  TC  18J 

164  IF  (kAI  l.NE.R  UlI  GO  TC  180 
RESKK  1=  RB  4L 
Ft S2  (K) *AP2H 
IF (K. EC. 2)  GJ  Tj  155 
K»2 

C  CHECkVir  nlfcV4tDUTC  LtrT  INTERSECT!  lM.tVALUATE  IF  RECLIPtO 
180  IF(lSL.NE.iO)  (j.)  T]  ibj 
I  Sl*l 
GO  TO  18 J 

160  If  ((F/JL-RBJL)/  (M1L-RGILI.CT  .0.01 
AA*(  AF2  Pj--aH  /  (PAJL-r.AILl 
«A-*A«cT-AV*:AlL 
Bb*  (AF2  IJ-.'F  »R)/  (KeJL-RHIL  I 
BB*  AP2.<~“  3*k  5l  L 
FESI  (K  I-  ir-r  A  1/  (  4  4—  T  B  ) 

RES2  (K)*AA*FeSl  (t.  J  «■  E  A 
IFlK.tJ.il  Jo  TJ  155 
K*2 


INTERSECTION, E VALuATE  IF  REQUIRED 


GO  TO  ItiJ 


INTO  <  1 10 
|NT  021*0 
INTC213G 
INT  02 1  40 
INT0215O 
INTC2U0 
I  NT  U2 1  TO 
INTC21EG 
JNTOilSO 
INTC2200 
IN TC 22  10 
IN  V Ji22U 
INT022j0 

Ini l224o 

INTC2250 
INTC2260 
IN TC22  JO 
INT022B0 
INTC22SC 
INT02JU0 
I N  T  02  J  10 
1NTC2J20 
INT02J JJ 
1NTC2J4C 
INT  02350 
INTC2J60 
INTC2J  JO 
IN  T  02 j8  J 
INTU2J9J 
INTC240J 
INT02410 
1 N  T  0  2  4  20 
INTC2420 
1NT0244J 
INTC2450 
IN  T024  tO 
INT024  70 
IN  T  C24.  fcO 
INTU249J 
INT025CJ 
INT  025 10 
INT 02 520 
INTC252J 
INT  02540 
INTC25  3 J 
IHTC25FC 
INT  02570 
INTC25B0 
INTC2  tSo 
1NT026>j0 
INUctlC 
INTO 26 20 
INTC26J0 
IN  Tl2C 40 
INT  026  50 
INT  02660 
IN  TC26 JO 
MTOOfctU 
INT  w26  VO 
INTC2  JlO 
JNT027 lo 
INTC2720 
INTC27JO 
I  NT  02  J  40 
1NTC2?50 
INT02760 
INTO  770 
INTC27EC 
INT 02  79o 
I N  TC2bOJ 
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FILE:  $«' VJft  F JRTRAN  PI  N  P  $ 

C  CHECk  F0«  DIRECT  HIT  U  -UOhT  JNTfejiiECTlCN.tVAUJ-rt  If  RECUIKEC 
163  IFUSY.Nt.lOl  G-J  Tl  161 
1SMI 
03  T)  It’ 

161  IFCRAl J.NL.R3IR1  CO  TO  U2 
RtSlI M  *Ril * 

RES21K1*  V*  <R 
IFlK.tw.ii  5:  TO  iaa 
K*2 

GO  TO  UJ 

C  CHECK  f  CP  i NT  tRMcCl AT  £  RIGHT  INT ER SECT  ION , EV U DATE  iF  REUUIRtD 

162  IF  1 1 SF.K E. 1C)  G  Z  TO  162 
ISR*  1 

60  TC  161 

162  IF  URAJ*~*  iJ’.) /<  P  UP-R31R).GT.0.01  uC  TO  163 
AA^A^FJ-ARZR)/  RAJ R-T4  1*  ) 

EA*AR2R~AA*n-tlR 

AB* (  A°  2  R  J-A7  |  /  (RHJP-R0. 10) 

EB  =  AR2  R-  AH  *  P  ?.  U 
RE  $UM*lAE-rtM  /  (  A  ■*— A  e  ) 

PES21K  )*  Aa*rcSUk  )♦** 
iF<K.Lc.2i  «C  TO  155 
K*i 

C  SE-1N1T  IAT  £ 

163  *»K2  h  J*«  i«2  F 
R  A  JL  *P  M  L 
FAjR=PAl» 

Pujc-Heii. 

R3jR*Pl  13 

1*0  C' IK'T  lf,'i| £ 

WKl It  16  .210) 

210  F3PH  at  (  •  LOG*  1>«0  THOJ'Jv»H,NO  PCINTS  FOUND* I 

6G  TC  la5 

IS*  IFIRCSII  1) .EC.  iCCC.c)  oO  TC  ISl 

If iresiui .e.>.svvC.j)  «>  to  isi 

If  IRES2U  )  .CN.SO.J  J.J  )  «)  TC  1*1 
1MRES2I2I  .fcO.iOuJ.U)  6  J  T3  IS  I 
GO  TC  155 

C  EVALLATICN  CF  1 KTERSECT1 CHS  when  CORVES  /RE  ALH1ST  TANGENT 
191  YGT»vj.0 

CAll  FcNKS  |N2,N1,X3,Y3,  23  ,tTl  ,YCO,  APRES  ) 

R1=ABSI  4Rhc3I2)-fpUS1U)  /2 

CALL  Ft'IFTi  (Ml,  I  t.A4.Y*.Zt  »HT  ^  .  TO  3.  A  RRE  S  ) 

«2*/eSl4RRtS12l--rHt<U)  )/2 

RES2(lJ=(  (-  i”*c)-Rc**2J*i  4*P2**2)  ) /12.0*ARP2) 

HES2 l2)*FrSZU » 

RE  ill  I)  =£  J«T(Ad)i  Rl**2>-«  RES2  11  )**2>)i 
RE§l(2)a-RcSl(l) 

155  RtTURN 
END 

SUBRCUT  l  It  PENOTjI  NA.'IY.A.  Y.  Z  ,  MT,  Y.J  .  X* ) 

01  PEAS  I CA  Y (4  0i . 2X <4 J).Z <40,40  ),Xl  *6  ),XPl lu> 

C  SEARCH*/ CK  J  Of  LOWER  Y  LINE 
CO  1C1  1*1. NY 
IF  1)111  »Gt . Y  j )  oJ  T  3  102 
10 1  CONTIMje 
1C2  JM-l 

C  CJAPSt  SCA I  FTP  ZERO  o\S3 

C  FE  XT  ♦  LlNcS  A*t  tXtOJTfcC  IN  FIRST  SEAR  CF  ONLY 
YC  * (  Yt—  YiJ)  >/i YtJFl)-Y(j) ) 

2X  U)*  in  1.  )•*  ll-Zll.J  )  )*Yr.+ZI  i,J) 

IFUNl  )  .  tC.nT )  u;  TO  107 
•I  0*2 

112  CC  1C3  I  =  I  C,  IX 

Zxm*(Zli  rjMl-ili  1  jl  )*  YC ♦ / < 1  .J) 

ZH*(ZU»  JUI-HI-  1.  J  ))*t  YO-YlJ))  /C  YI  JFD-YI  J) )  F21I-1  ,J) 

If  (Zxt  1  )  .LC.hT)  uC  TC  US 

IF  UZXH-D-HT)  /Izmi-hT)  .LT.O.J)  GC  TO  U  , 

103  CONTINUE 


INTCZulO 
INT02620 
IN  T  C  2  S  20 
I NT  02  6*0 
I  NTOJH50 
1ST  02  t’c.0 
IN  To  Co  /J 
iNTCZtfC 
In  rc^tisu 

I  NT  02 SCO 
1NTC25  1U 
INTC2S2J 
IMT  02‘>  jJ 
IN  T  C2S4C 
IN T 02S 50 
1 1 J  Tv2  SoO 
1NTC2S70 
|NTw2SC0 
IN  T £2 6  70 
INTOlOoJ 
INTOjOIO 
INTC3C20 

In  t 0  jo  j o 

INT 03040 
INTL3C50 
INT  030t0 
IN  10^070 
INTC  30Eo 
INTDOovO 
IHT03100 
IMT03  1  10 
1NT  J3l2o 
I N  TC  J I  50 
INTO  21 40 
INTO 31 50 
iMTOiltO 
INTO  11  70 
INTC31SJ 
INTO 3  ISO 
1NTJ320J 
IN  T  C 12  10 
IN  TO 3220 
IN  7012 10 
IN  V  C 12  40 
1 NT 032 SO 
IN  T  C32 1 C 
INTO 12  70 
INTC126J 
INTC22S0 
INT013C0 
INTO  3l U 
INTC 2220 
IN  TO! 3 30 
1NTC33H0 
IN  T C 3 350 
I  NT  O31t<0 
INT  C 13  70 
INT 01260 
INT 03340 
IN  T C  3400 
INTC14  10 
INTO  14  20 
INTC143C 
IN  T034h0 
I  111  0.1450 
INT  C34CU 
1  NT  C04  70 
INTL34EC 
INTJ1490 
IN  TO  j  5  00 
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FILE;  SWVOR 


FORTRAN  PI 


N 


GO  TG  lit 

C  CALCULAT  iC'J  2f  PFM  ET?  Af  I  f)N$  I*  E  I  FFC*  TRIANGLE 

104  XK*xu-nMxu)-xu-m*irc-v<jiJ/mJ*ii-T<Jii  • 

IF I2-.EJ.H  n  GJ  n  IC5  .  _ 

IF  1(2X11 1-nT  )  /  t  Z'-'-RT  J  »LT  .0  .j )  GOTO  1C5 
00*AE5th  WV)  iaoS  (ZX  tl-U-HT  1 
X$*>t  l-ll*tX'i-Xt  l-ll  ZXt I-1I-HTI/0D 

CO  TC  110 

ICS  >S«>PMXUJ-XM)*AdSUP  -hf  J /( AbS (  ZM-hT  i  ♦  AES  t  ZX  ( 1 1-  FT  1 1 
bC  tC  110 

C  LCbl  Cal  FLI\tl*ATiCNS  ACCLRULAT 1CN 


10  7 

>R  ( 1 1*  X  (  1) 

1*1 

GO  TO  113 

loe 

XR 1 P  )*X( 11 

GO  TO  111 

11c 

XR  ( 1  R I  *  X  S 

111 

1FUC.E0.2)  GO 

113 

1K*2 

16*1*1 

GC  TO  *12 

116 

IFIiR.tC.il  CO 
XR(2)*1UC.  0 

GO  TO  114 

117 

>8111*1000.0 

XR  (  21  *  l  CCC.  C 

114 

RETURN 

ENL 

SURRCL  I  INF  HTPLTU,V,Z,XifY,;,M,NYnPcS) 

C  LINEAR  IKTcSFClAMCN  CM  C  AL  i  F-  AT  l  IN  SU- FACE  TO  EVALUATE 
C  CP  VALLE  t  T  CO- )r\Ol.N«  ftS  XGtY.,.  TmE  RESULT  RETURItO  IS  ZFES 

ni^fcAMCN  A(<tu).  VUii)tI(NO«40) 

C  SEARCH  FCh  1  CF  Ltff  X  L  INt 
JO  1  1C  *  if  N  X 

IF  (X (  IC )  .Ll.Xu)  GO  TO  2 

1  CCMIME 

2  1*101 

C  SEARCH  FCP  J  -if  lower  y  line 
oo  j  o;«i,i\y 
IFCYtJC J.SE.fG)  GO  TC  4 

3  CCFT1FJL 

4  J=JC-l 

ii*i»u*i)-Mji)/uu«i)-Atm 
BL*Y t  Ji —At*  X  I  1 J 
YC R*  (  AL*  Xj  i  »8L 
IF  (Y C  .  EO  .Y  lR  »  03  T1  1 

iMtj.ul.y.K)  jj  t:  5 

C  RESULT  1*1  G  Z  IS  I  I  LUwRtK  TRIANGLE 
X3 *X  tin  I 
Y3*Y( J) 

23*1 ( It  if J I 
bC  TC  6 

C  RASLLTINu  l  IS  IN  UPPER  TRIANGLE 

5  X3*Xtl) 

»3*Y(J*l J 

C  CENE^ALZC#Ic0LAT ICN: CCCD  FOR  EOT  F  TRIANGLES 

t  >i*xtn 

Yl«Yt4I 
21  *2  1 1  t  J  J 
X2*  Xl  I ►  1 1 
V2*Y (J ♦  l » 

-AP*t ll ;  1-2  if*/ YI-V2I-UI-Z2IM YI-YjI I/ICXI-X3JMYI-Y2 J-UI-X2 I 
1*CVI-YJ i I 

tP» (  (Z l-Z2)-4P*lXl-X2I»/CVi-Y2 I 
CP*  Z  1- 4R  *  X  1— •'P  *  Y  1 
2KES*AF*XLYrtF«rw«C’ 

GO  IJ  100 

C  RESULT  INC-  2  IS  O  t  OIVUING  LINE 


iNTCJSIC 
N TO  3520 
NTCJ5JC 

{NT03540 
NTCJ550 
NTCJScO 

iNT035  10 
NTCJ5E0 
NTC3590 
NT03600 
IN  1C  36  10 
INTO  3020 
INTO jbjO 
INTJ3o40 
iNT  J ib50 
lMT03t>60 
INT  03c  70 
INTOAbbO 
INTC36S0 
IN  TuE/JJ 
I  N  T  0  J  7 1 J 
1NTC3720 
lNT037aU 
I N  T  0  3  2  40 
INTC3750 
INT  03  1  tJ 
IN  T0J770 
INT03780 
INT03/9J 
INTC3ECC 
II  TO  Jd  10 
1 NTCJB20 
IN TC 3630 
INT 03 843 
INTO.)  tl  50 
INTC3660 
INT  03  8  70 
1 M T  ujb 3  0 

INTC3SS0 
IKTU390J 
INT C  3  5 1 0 
INI  03320 
1 N  TC j  93 J 
IN  TL  3540 
JNT0395J 
INT0b960 
NT  039 70 
NT  039  60 
NTC3550 
NT  C**olO 
INTO401J 
INTC4C2C 
INT0402J 
1NT04040 
MTCtCSJ 
,  NT 040o J 
INTC4C7G 
INT04U90 
INTJ4O50 
INTC41C0 
INTC4I1J 
IN  T04 1  ^0 
IMTC413C 
I.nto^iaj 
INT04laJ 
IN  VC41L0 
INT0417J 
INTC41EC 
INT0419J 
I  NT  042  00 
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APPENDIX  3 


LISTING  OF  WVDR 


j 


] 

\ 

i 

j 

i 


U/Lt/i*  If  .4 


S3 


Flits  tVLH  F  CRT  KAN  FI  N  P  S 


C  MARI  OF  HUT 

JN7000  1.) 

Cl  V£  NS  IN  ilF  (1J),’H1(  lol  ,  *\(  Hl.yil  4C»  , 

INI  T  J 0020 

1  XA  IT  )  ,  ,0  J  ,<  3(40  1  ,  Y  J  IaO  ),  2F  t-*0,  4u  )  ,X  AXt  AO  J ,  Y\X1  40)  ,• 

II  TCcCiO 

2Z  A  A  l  ^  L  #  U)  t  A  hO  t  V  i  Jt  l  H  0)  fi  viX  140  t^o) 

mjj.no 

119 

tCRU T  1  ,  , 

INTO.C5C 

131 

n'^MiF/.n 

IN i OoOoO 

no 

f P IAT( t9.ll 

I  NT  000 7 J 

120 

INTCOJfco 

132 

fj  <«Aur  ic.  si 

IN  i  C  C  C  S  w 

135 

F  ,)r- 'NT  1  K  la  .5  1 

INTOolJO 

13  < 

F!>N  *r  I»-l  /  .5  ) 

INTlCLIU 

C  READ  ”0  IFF  '■.linarijN  U  AT  A 

i'lT  00120 

R  C  A  C  (2i  119  1  U,  MY 

MTJOliJ 

C'i  121  l»i  (i>< 

I‘J  1  CO  1  -,0 

Rfc'.UU,  liJi  AA(J) 

I Ul o  J i jO 

121 

CC NT  IN  IL 

IMTCCUC 

01  1  tt  1*1,0 

INT0J1  /O 

REJC(2,nJ<  Y  A  (  1  1 

IMTOol A j 

122 

CC.ST  IN  Jt 

I  *4  T  L  W  1  ^  J 

0  1  12  J  1=1 »NX 

IN  T  002  00 

CJ  123  J  =  i ,  l» 

IMco2lO 

R  E  A  0  4  2 ,12  01  .At  1  ,j) 

111  f  0  J  £  .  J 

123 

CUlMcE 

INTO  02  30 

CJ  221  1=  i,  IX 

I  NT:j2-»J 

REaDI 2,130)  XJtlJ 

INI  C  u  2 : 0 

221 

C.TI  IT  IN'Jg 

In i  jw>2oj 

OC  222  i  =  1 , NY 

INTCC2  70 

RE  AO  t  2 , 1 i u)  Y3I1J 

I ‘IT  J02  3J 

222 

CCNT  IN  IE 

IN  TO 02 00 

OU  223  1*1  , NX 

IN  T  l  C  3  CO 

0(1  223  J=l,'l» 

INTOJoU 

PEAL (2 , 1 2 J i  Z  31  1 , J  ) 

1 II  T 00320 

22  3 

CCMlIMt 

IN  TO J 3 30 

C  TRAMSSTA^  '.iin1,  A  riJN  1A1RICHS 

1NT003 a J 

CO  224  1*1  ,N.Y 

INICOiiU 

224 

*\A( l)=Y4(t | 

iNTCOitO 

CO  2.5  1=  1,  'x 

IMTCJ3  lo 

22  S 

YAM  i  1  *Xi  U  I 

INTlCicG 

DO  iic.  I  *1  ,  JY 

INI  00390 

UO  226  j  =  i. .  I< 

1NTCO400 

221 

2A>( 1  , w ) * 2 a ( w  » 1 J 

i'ITU.410 

00  2.7  1  =  1,  IY 

I  NT  00420 

22  r 

XbX 1 1  I  *  Y  E  1 1 1 

INTI  1.4  30 

00  i  eb  1=1,  NX 

IN  f 00a 40 

22  a 

YRX(  1  J  =  Xfc  m 

[NTCJ450 

CU  2  29  1  =  1  ,NY 

INTCC460 

DJ  2<9  J*l, Ia 

I*lTuU4  70 

229 

2oX  (  1,  J  1  *//(.!,  n 

I N  T  004 oJ 

C  FFiCBt  StniNuS 

IMTCC490 

YP  A*  0  *0 

IN  T  005  aO 

YPL<=55. 

IN  7005  10 

XPA=C.O 

INTC0S20 

XPC-2C  . 

IN|TJ05J0 

XfO=-2 J. 

INTCC54C 

C  PR  CG  RAF  CO.vSIViTS 

IN  TOO 5 50 

YNl.s— jJ • 

IriTCOOjO 

YREF=30. 

I’l  tCa5  70 

YR  1  '4=-3  0. 

IN.  TOO  5  00 

XPlO-oO  . 

*  INTOaSRJ 

>RIP=3  0. 

1  Ml  COc 00 

XR  IN  =—  3  o  • 

I  NT  JviO  ll 

ICCFS=i 

1  NT  0.6  20 

.NJ0upS=2 

INTOCtSO 

NO  S  1  '•'=  1 

1MT0J6AJ 

i  i  r  *  i 

IN  TCCl 5C 

I  SC  AN'=  1 

IN  T Goc  c>0 

R£L  X  FS=  J  .  J 

IN  TOON.  70 

C  EXFfcFl NcNT  Si'LLATiCN 

HTCctEC 

19o 

C".htp'u: 

INTOoSNO 

401 

ftrNAT  (  r  15  l 

1NTCC  7C0 

F  HE: 


FORTRAN  P  I  • 


402 


403 
C 


50  C 


READ! 5,4011  PT 
WHITE  (o.-*uli  ‘•T 
If  |NCM  M.  tM.2 I  tC  TC  500 
READ! 5, hC1 I  >>s 
wtft  1 T  i  loi'tul  J  °5 
St  AO  (5  *hCi  I  fLI  A 
*R  I T r I  a, 401)  ALF ». 

FEACra.YJli  RFii 
SR  life  I  a  .'lOll  PrU  i 

CALL  i  I  r  PL  i  (aA,  r  A  ,Z  A  ,ALFA  ,Prtl  I  ,NX,NV,CPAI 
AlFC=AlFA-XPC 

CALL  MrRLlIXA.rA.ZA.ALFC.PHII  ,\X,NY,CPC) 
ALFC=ALEA-XPJ 

Call  i  i\  t  r  l  i  i  x  a  ,  y  a  »  z  a  ,  a lf  o  *  pm i  i  »  nx  »  *-iy  ,  cpo  l 

PHIP=PH1 1-YP  - 

CALL  1MT  °LT  ( •<  ° *  Y  c  1 1  r- .  Mr  4  ,  =H  IS  ,N  A,  ‘I  Y  ,  CP  A  ) 
WRlTEli>»VJt)  CP^tCrtiCPCfLflJ 
FORMAT!  •  T-IJL  CPi  *.  4r  10.  5) 

FoYf<  -  FT  -PS 
PPA*CFA*PO YN+PS 
PP4*C Pe*°C»  '+pi 
FFC  =  CFw*f  Jff.tPS 
PPU  =  CP J*PaY  l*-PS 
WR I  T  E  14»  »4  1 1 1  PFA,FFE,PPC,?PD 
FJ-!MAH»  ii'KL  RES'  »4F 10. ?  I 
RE  PC  "bPSUAt'M  ‘T  3  CAVA 
IP IXCSiV.l  ,.l  I  oC 


*0  501 


181 

TO  30  5 


REACH  i,  12£i  ••p. 

W  fi  IT  E  ( a  ,  I’-'i 

PEAC15.13CJ  PPo 

wRireio,  lioj  pps 

Rk<Cli*iJ')|  PPC 
WR I TE ( o  »  1  2 c I  Pr'C 
RESClb.ljai  -’PJ 
WRITE  16  .lloj  P»0 
REA  J  (it  13  ;  i  Pii.J 
«RlTf(o,ij6J  P-.1N 
lFIISCAHLO.il  aO 
IF  (  Ht'A'l  .01  .  ICO-  I 
1 1  T*  1 
ICCPS*l 
ISCAN=  ISCV-*  1 
FSl  MF31H5.  J 
WRITE (t,  13 tJ  Pi  10 
CCNTlf'JE 
F$=PJ1A 
P5C=PS 
FST*  PS 

C  CALCLLULS  PRESSURE  COEFFICIENTS 
300  COnTiriF 

WR  1  T  t  4  j  r'f  'iij  I  l  IT 
FuRMATI  'C  ,*  1  It -AfiCN  HO.  • 
PUY  l=PT -Pi 
CPA=( PrA-r  >) /°0Y  * 

CPJ=(Pp  I- >i)/P jt  . 

C FC=  IPPC-P.R  I/’  Of  I 
CPO-l  •VJ-FS)  / r.i Y i. 
^!iCIj,/U  cP4,c.>-nr.pc»CPii 
912  FOM’"lf  oLCJ  Ll-S*,4flJ.f>i 
C  CORRECTS  P  j  -ijiL^PTHi't  TO  E  i  jjr.c 
IF  1NCC  iPj.tl.tl  il  T1  aJl 
2AFAX*- l £000. 

.CALL  "ur.>kiftr',,*i  IX  ,-VY  ) 
IF  <Cn  J  .CT  .0  1  14  X  l  p‘,0-  1  .0  i*i-3 
If  IC  PC.  J  1.  L  -  /..<)  r  :  0  =  1  •  )  l  X  F S 
IFICPU.Ol  .C  i  'A X  1  p  jC  =  1.  j  1'PS 
ZfcVAX  =  -l  J./ j  J  • 

CALL  vxA'lrf  li->)lt»<Xt\X.>:Y) 
IF IC°E .«T  ./ 2 1 A X  I  PjC^i.Oi^Pi 
ZAPl MIOOJ  ). 


£01 
18  C 


181 


94C 


908 


.151 


CP  ARE  IN  C  ALI  3R  AT  I'JN  R  ANOE 


ri  r  jo  /  10 

INTCJ720 
IN  T  CO  7  20 
IN T 00  7  40 
IN  1 0  j  7  50 
IN  T  C  U  7 1.0 
INT  Jo 7  70 
IN T C 0  I E J 
INT to /So 
I N  TO  J  9  00 
IN  T  L  C 1 1 0 
IN  YOOc*  20 
IN  T  0  0  L  j  J 
INI  TCoc 4 J 
INToCiPJ 

I N t  uuooO 

INT0o6  70 
1  NT  0  08  80 
IN  TCtcSO 
I NT o  0  5  J  0 
INTO  JO  lO 
IN  TLCS2C 
HT00N3J 
I li LvSXO 
lNT0o950 
IN Tcoo 60 
1NTCCS70 
lNTUONdO 
1  N  TO 09 90 
I ■'!  i  C  loCO 
INYO  10  10 
I  iTUCZG 
1NTG10  20 
1  i| TCI 0  .0 


I N  T  C  1C  £0 
INTO  1060 
INTOIO  70 
INTO  1080 
INTO  10  NO 
IN  TO  11  CO 
I’lTClilO 
INTOil  A» 
INTCllJO 
1NT0114J 
IN TOll 80 
INTCll/C 
INTOll/j 
I'lTOliaJ 
INTO  l  1  *o 

Ini  On. oo 
IN  T  C  12  10 
1NTv122o 
INT  U230 
IN  TCltiO 
INT  Jl2  50 
I N TO  1 8a J 
nrci2  jo 

INTOitoO 
I  I  T C  1850 
INTO  i.oo 
I  - 1  T  0 1  a  1 J 
IN  T  o  1 J  20 
INTO  1  i  JJ 
1  NT  0  la HO 

IN  ToliiJ 
I  NT  013  6 J 
IN  To  l  2  JO 
IN  Y  0 1  2  fc  J 
1  NTOi  j'Jo 
IMTCHCO 


73 


FILE :  WVCF 


FORTRAN  PI 


CALL  HI  (ARYlZA.ZAHIN  ,l\Xf  AY) 
IFtCPA.Ll.ZAIlJ)  PsC»«.99*PS 
IF  (CPC.  LT.  t  A*  (  \| )  F9C*0.99*PS 
_l  ic  «C.99*  PS 


963 

966 


S2C 

921 


.13  i 
TC  ISO 


' »  F20 • 5  ) 


IFICFU.L  i.ZAItiU  P5C 
ze^iN*  u jjj  . 

CALL  PiNAFY(2e»ZE**NtNX,KVI 
1  Ft  CP"  ,L  r  .2JM1N)  P9C»C.99*PS 
1FIPSC. EC. P3T)  v»C  TC  3U1 
IC  CPS*1 C  C  F  S* l 
*RlTEto,9fc4l  IC1PS 
FCFMTP  1CCPS  i9* 

IF  ( i  CCPS.OT. 1C)  oC 
FS*P$C 
FST*F$C 

WR llEto. S21)  PSC 
FjRMATP  °  S  C  IS 
OU  ID  300 
301  CCKTU’UE 

C  CALCLlfTtS  t>.D  ALF  f.  Pnl  A  N  jL  fc  S  FOP  'I*  £  'll' 

CALL  IN  UCSt  VP  a.  yP3  ,VRl.J  ,YRUP,  YRI  N, 

1CPA.CP2, Al F,PHi, *4, YA.ZA, 

2>B  .YrUzd.XAT,  Y4X  ,xAX  ,X  EX  ,Y  EX  ,Z  eX.NX  ,FY  ) 

AL  r  1=AlF  t  l ) 

ALF2=ALF(2) 

PH  I  I  *Prl  1(1) 

<>Hl2*P8i  (2) 

C  CALCULATES  TwC  4lF  £  Fn i  ANCLES  FIR  'I'  £  'III* 

CALL  I  1 T  SCSI  APA  »XPC  ,  >RLQ,xr«UP,X*IN, 

1CFA, CPC, Pn It  'L F i a 4X i Y AX ,Z A  X, 
2>A»,YAX,Z4X,XA,Y*,ZA.X4,YA,ZA,FY,NX) 

ALF3«ALf  ( 1) 

ALF6-ALK2i 
PH  1  3  *  Prl  I  »  1) 

FH14=PHI  U) 

C  SELECTS  Tnc  FF  .iPcS  ALF  £  Pnl  INGLES  HUT  OF  THE  FOUR  VALUES 
HALF* I  (  ALFi»  \LF  2  ♦  ALH  3* -»LF  4  )- AM 4X1  t  A  LF 1  ,  A  LF2  ,  ALF3  , 

1  ALF4 )“A '  I  !i  (4i.fi  .AuF2t-*LFwfALF  6)1/2 
FPHI*(  (Pnil  ♦Fni^*Pni4+P:iI<.)-AM4Xlt  Prill,PhI2,PHl3,PF(6) 
1-AM  IN  It  P  F  I  l.PHlc.Pill  i  ,  PH  I  4)  )  /2 
900  IfHITEto  ,)0i)FALF,FPhi 

SCI  F  JRPATt  '  SLltCTtO*  ,2F10.2) 

C  CALCULATES  lEw  Pj  FR'H  DATA  OF  POSITION  D  A  NO  CORRECTS  PT 
ALfCMraLF-XPO 

CALL  lNTFLTtXA,f4,ZA , A LFO A ,f Phi , NX , NY »C PO) 
ALFCMFoLF-X’C 

CALL  1  \  T  F  L  T ( X  A  •  Y  A  .ZA.ALFC  1  ,F PH  I. NX , MY . C PCI 
90S  n'RITFtofVlJ)  CPC, CPU 

910  FC  F*  AT ( '  JEJ  C’C  £  C°C  CALLEJ  «.F 15. 5) 
PSN«IPi>C*CRJ-FPJ*CPC»/  (LFJ-CFL) 

907  WUTKb,,UI  FAu)  .nPHI  ,P  3  t  ,>T 
311  FCFPAT(*  ITi-iAT  (CN  VALUES’ .2F6.2.2F10. 2) 


n 

m 

S5i 

S6C 

961 


O.NVfcxi.EME  iEoTS 
CtLPS=PS  A-RS 


HKlTct  j.9'32)  PS.FS'I  fCtuPS 
FORM 41(«  tPSRS  VALUE  S' »3F 15. SI 
E  p  s  p  s*  a  & ;  ( t»vi-js  i/p  s  j 


FKlino.itl)  EPjfS 
FORMAT  (  •  cR  S*>  S  !„•  .  r  25.  15) 

tPSPSo*.  )JVJJ JJlh.099 '(PT**2 .o  )- .VJ00009B  *RT 
IFt  IcFjl.EO.  2)  c  R $PSG*.  L COOCC  JC\)09*  t  FT**2.J)-. JU009Ud*PT 
IF  t  ltPj<  . E  ) .2 )  tP  jPSiia.COJu^bV(.Oj9'(  PT"Z,  JJ  — .  C  w6vJ  w  C  6*  P  T 
IF  IEFSPS.uT.EPSPSl)  EC  TO  310 
-  CO  TO  i j * 

CCFTIFJE 
II(*IIlH 
IF  II  IT.  GT.  10)  «0  TO  305 
FS»PS*RtLXFS*(PSA-PS) 

PSOOJ, 

PST-PS 

S2A  Ml  It  16, 925)  PS 


310 


IN  To  14  10 
INT0192J 
I  N  T£ 1420 
IN  f  0 W40 
INT 0145) 
1NTC146C 
I NTO 14  70 
INTC14E0 
INT  C 1490 
INT015JJ 
INTC151G 
INT0la20 
lNTOi.539 
IN  TC1540 
INTO  15  5.) 
INTC156J 
INTO  15  70 
INTJ15SJ 
I  N  Tt  1590 
INTO  loOO 
INT01610 
IN  TC  1620 
INT 01 630 
I N  Tu  1640 
INT  C 1  £  50 
I  NT 0ifc60 
INTClfc'iC 
IN  T  0  lb  00 
1 N  TCI 690 
INTO  1  l  CO 
INT017 10 
INTtl 72 0 
IN  To  1 720 
I  N  T01 740 
IN  TC 1 750 
INTO l 760 
I  N  TO  1 770 
INT  C 17«0 
INTO  1790 
INTCieCO 
INTO  18  10 
INTJlt)  20 
IN  T  Cl  830 
INTO  1840 
INT01B5U 
I JTCietO 
INTO  18 70 
INTClbdJ 
IN  TO  1890 
INT01900 
INTCIS10 
INT Cl  920 
INTO  19  30 
INTCIS4C 
INTO  19 50 
I N  TC 1930 
INTO  19 70 
I  NT  019 00 
IN  TCI 990 
INT020CO 
INTwZOlJ 
INTC2C20 
INT 029 30 
IN  TC2C-.0 
INT 020 59 
1  NT 02 060 
1NTC2C7C 
IN  TO 20  <)0 
I NT02090 
INT02100 
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Flits  SVC? 


FCRTRAN  PI 


25  FORMAT! 1  CURK  PS  IS  *^15.51 
GO  *0  300 

03  J  f  ALF,FPFi*P3'l,PT 

71  FORMAT!  JF  la.  2,F  15.  J) 

cj  re  in- 

C5  UR1  TE  (r>  ,107) 

07  Fi)?-*AT!»  iCA'I/EJ  100  IIMES'I 

GO  T0  WJ 
E W 

SURROUTI'Jc  l  iTSC:>|:.RP1,ARP2,.lU0.4RUP,ARIN, 

1HT 1  ,HT 2  j«c3 it  PcSitX4»y !•  At  » 

2*2.Yx.2<*X 3tY3.2J,x*»,Y4,Z4.Nl  ,A2I  ,,  ,  „ , , 

U1 MfcAs  i\N  RlS1!IuJ,=  ES2I  10),^.nkc51  10), *1140)  ,Yll40l  t 
111 14C,4C) ,*21401 ,V2 1401,22  t4J  ,40  I  ,X  J  140  )  ,Y  j!  40  1 , 

223! 4C»  «  C  )  * *41 40) ,  V ♦! hC) *24 1 40 , 40) 

PtSi (1 )  = too  J. 

PE  II 12)  s^O-jO. 

RES21 1 7* ijJo. 

ReS2  12) =9000. 

ISl*l 
1SR*  1 
0AR=1. 

AR2R-ARIN 
GO  TC  15* 

5C  A«2ft*A.-(2K*dAR 
AR2RJ=AR*R 
52  ARH*AS2F-ASF1 
ARM2*Ak2R- W2 

CHECK  FCF  LC«fc  C«LIERaTIO!  RANGE 
IF  l/Pwl.lT«ARLO)  GO  TC  150 
IFIAPM2.LT. A^lJ)  GJ  TO  150 
K*l 

CALCLLATtI  INITIAL  PcNGTFATICN  PCI  MS  , 

CALL  PtMPlM  11,14.;  ,*1  ,7  1,l  1  ,HT1 ,  ARMl  ,  ARRESi 
PAJL*A.»Pti  ill 

RA  JR*A 4 RE  S i 2 )  .... 

CALL  Pt'IPT  j  II  1,I2,X2,Y2,Z2,HT2,AR42,ARRES) 

R8  «l =Ar  FLS  11  ) 

RBJR=AA?.ESl  21 

CALCULATES  SUCC  £3  1VC  PEN  ETRAT  ION  «>C1NTS 
DO  140  *-1,10000 
AR*R«AUROAR 
ARMi=AA2R-lA>Ji 
ARM2^AP2R-A<32 

CHECK  FOP  ijP>ER  CALIBRATION  RANjE 
IF  !APf‘L.uT  .  A~  UP  I  0)  TC  195 

IFIARP2.0T.mkUP)  uO  TC  19a  . 

CALL  Pt’IPT  il'U,  U,X1,  tl,2  1  *  H  T  l  ,  ARM1  »ARRE  S) 

RAjl=ASFES 111 

CALL*  PC 'IP  fS  !  I  l,‘12,X2«Y2,Z  2,HT2 , AR 42 , ARRE  S) 

FBI L*7SScS  ti I 
RBIR=»ASRES!  21 
1F1RAJL .EO.lddJ. J  ISfiO 
IF  IRfcJL.t  ..1000. J  1 3  l  *  10 
IF  CR  PJR  .t) . 1 JJJ.  1  1SR*10 
IF  1SA0R  •  EG  .  iOOd  •  I  I>S=lu 
IFIRmIl.E.. 1JCL. I  131=10 
IF1RPIL.EO.1jOJ.  i  I3,_=ld 

If iPeiP.Ec.io >0. i  iss=io 

IF  l9  AIR  .  EO.  1000.  I  i  J':  =  ld 

CHECK  FCP  CI4ECT  FIT  J1  LEFT  INTER  SECTION, EVALUATE  IF  REOUIRED 
IFtllL.Nt.10)  Ou  TO  lo4 
-1SL*  i 
C-C  TC  180 

E4  lFirAlL.Nt.ROlL)  GO  TC  1B0 
RE  SI  IK  )=R‘*  IL 
SES2  1M  *7r2S 
IF1K.E..2)  G.  TJ  155 
K*  2 


iv  r  o  <  i  io 

I MTOxl 20 
IS TL212C 
1 JV02140 
INT0* 1 50 
INTC2U0 
IMT02170 
INTOilEG 
INTC21S0 
INT0220J 
1NTC2210 
IN  T 022*0 
1 NTU22  o J 
IN  TC22  40 
1MT0225J 
INT  022o0 
INT  022  70 
I  NT  022  00 
INTC22S3 
INTC23G0 
IVT02J10 
IN  TC2320 
INTC2.>30 

111  TG234 J 

INTC2350 
INT02O60 
INTC237C 
IN  i  j  2  3  bd 
I  N  T  C2 .»  90 
IN  TC 24 CO 
1NT02410 
1NT02420 
INTC24iJ 
INTu2440 
1NTC2453 
IN  Y  C2m  c  J 
INT 024  70 
1NTC24EC 
1NT02m‘5J 
I  NT 025  JO 
IN  T 02  5  10 
IMT02520 
1NT02530 
IN 10*540 
INT0255J 
INT  025. >J 
INTC2570 
INT025B0 
INTC25SC 
IN  V02tOO 
lNT02t,10 
IN  T 02(  20 
INT02O0J 
INT02G40 
INT02050 
l NT  02060 
INTC2C  70 
INTC21E0 
INT02O90 
INTC27CO 
INT 02/10 
INTC2720 
1NTC2  7  30 
INT  127  40 
INTC2750 
IN To 2 /fed 
I  NT  02  7  7  J 
IN  TC2 7tO 
INTC27S0 
INTJ2800 
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FILE:  WVOR 


FORTRAN  PI 


N 


161 


Gil  TO  Ifl j  ,,,M 

C  CHECK  ICR  itfERMEUIATE  LEFT  INTERSECTION, EVALUATE  IF  REQLIREJ 
180  IF  USL.Nt.lO)  GO  TJ  16U 
ISL»  1 
CO  TO  i 33 

16  C  IF  (I  F>.JL-Rfc3L)/l  1  AiL-SrilL)  .GT.O.OI  GCT1  Id  3 
AA*  t  AF  2R  J-  ‘A  HR  I  /  (  K4jL-r.A  I L I 
E A*AP* R- / A$R  AIL 
AB  =  (  AK2  3  J-a  «2k)  /  iPdJL-PslL) 

BB*AR2R-AU.*FRIL 
KESl  IK) =|db-FA)/  l/J-<ei 
RE j2 Ik) =AA^RLil( K) ♦HA 
IF  IK.LC.,2)  G  1  l"1  Ip5> 

K*2 

C  CHECK  F  JR  DIRECT  HIT  J1  SIGHT  I hTE RSEC TI CK , E VALUAT t  IF  F  ECU  I RED 
183  IFI1SF  .Kfc.iO)  GO  TO  loi 
I  Sk*1 
GO  TO  163 

IF  IHAl4.KE.ReiR)  CO  TC  182 
KEilIK)=i'H  A 
PES2  ik  .  =  axh  r 
IFIK.EO. 21  uJ  TJ  155 
K*  2 

GO  TO  16  3 

C  Check  FOR  IinTEkMEJI  ATc  FluriT  intersection, evaluate  if  reojired 
182  1FIISF.NE.1J)  G3  TO  lt2 
1  SR*i 
GO  TO  Ui 

IF  ((RAJ  F-RdJ  •'.I/IRAIR-RBIR )  .GT . C.O )  GC  TO  163 
AA*IAF2R--A««.f<J/IFJJH-HAlR| 

BA-*  AR2R- AA  AIR 

AB=  (AF23.J-AR2R  )/  I  5  3J  R— SB  II) 

8B= AR2  R-A3+RilR 
RESi(K)=IHA-HN)/I AA-AB) 

RfcS2  IK)=AA*PcSilK)»£A 
IFIK.FJ.2)  uO  TJ  155 
K*  2 

C  RE-1M  TI  ATE 
163  ARiRJ  =  AR2R 
RAJL  =  F  AIL 
R  A  „R 

RB  JL  =  P;3  1L 
Re«K*Fblr 
CONTI  Nlc 
HR  ITEIu,  21  )) 

FGFKJT1'  L'JOP  140  T  FPGU  On,  "10  POINTS  FOUND') 

GO  10  155 

IF  IRES*  II  )  .E0.9JUJ.U  I  GOTO  IV  l 


162 


1 A  C 
210 


195 


IF  tHtSl  12)  .e). 9000.0 1  GO 


191 


IFIRFS2I  1)  .l  J.  KuJ.  0)  GJ  TJ  191 
IFl»fcS2  121  .t  J.VOGJ.d)  uJ  TO  191 
GO  TC  155 

C  EVALUATICA  IF  INTER6E  Cfl  u‘IS  WHEN  CLRVES  ARE  ALMOST  TANGENT 

191  YGC=0.J 

call  PCNPTS1 N2,.Nl ,  X  3  ,  Y  3  .  Z 3  ,H  TI  ,YG3,APRES) 

R1  =  1  tS  I  A  F  r.  l  S  I H )  -  A  R  ;  £  5  (  j.jj/2 

CALL  ♦  t  K  F  T  >  I  v2  ,n  1  ,  <4  , Y4  » i  ♦  »nT2, YuJ,AFRES  ) 

R2* A’St  ARR  f  SI  2  )  —  •-•  *'  R  l  51  1)  )  /2 

FES2  II  )=  (  l*-  1**2  ,-(  <..**4  )♦(  A :<n  2 ^ )  )/I  2.  0*AAP«:) 

RE  J2I2)  =RL>21  1) 

RES  1  (  1  )=S<)-.TIAio(  (31**2)-|><tS2l  l)  **  2  )  I  ) 

FESl  12)  =  -  A  t.3  i  1 1  1 
15!  .RETLsN 
ENL 

SUllF  CUT  I  Nt  F  CNPT  S I  NX , NY  *X  ,Y  »  Z  »  FT  »  Y  G»  X  R  ) 

C 1  ** £ N S  I  )N  Yl  4C)  »lXI  40  ,Z{  40,40)  ,AI40)  ,XPI10) 

JR*l 

C  SEARCH  FCF  J  LF  lCKEc  y  line 
CIO  lul  1=1,  NY 
IF  I V I  ll.Lc.Yu)  uC  TO  10^ 


INTUcJ  10 

I.NTC26  JO 
INTJ^j^j 
I N  T  0  2  b  3d 
IMTi.JFf.'J 
INToio/J 
INTCZeeO 
INT 02140 
INT02500 

in  re.uo 

IN  T  02920 
I  NT 02 9 30 
IN  TC2S40 
INTu295J 
IMTCHStO 
INTJ2v  Ju 
I  NT 026  30 
IN  TC 24 SO 
INTOjOlO 
INT 0  301 J 
INT03C<0 
A NT u 3030 
I N  TO j04 0 
IN  T  C  30  50 
INT0J060 
INTL3C  ) 0 
INT  JiuHG 
1 N  T  CO  090 
IN  T 1 3 1  CO 
INT031 U 
1NTC512J 
IN  TO  3  1 30 
INT0314J 
I  N  TC  31  50 
INT0316J 
I  N  T 0  3 1  70 
INTC31EU 
IN  T 63 1  90 
I  N  T  0  j2  C  J 
INT  C  1 0 
INTO  1220 
I  f •  T  C  l  2  3  0 
IN  TO  3  2  40 
I  N  T  6  J  2  5  0 
IN  TC J2t 0 
I N  T  u  3  2  7  J 
I  il  T  0  3  2  30 
INTu j290 
INTO  33 00 
INTC33  10 
INTO J32U 
1NTJ  i  i JO 
I N  T  C  3  3  *»  J 
INT  03350 
I N  T  03  3  6  0 
IN  TC  J3  70 
INT0333U 
INT0334J 
INT00400 
1NT03410 
INTC3420 
INTCj430 
I  N  TO  j4 40 


IN  TC343C 
INTO j4cd 
IN  TO 34  70 
INTC34E0 
IN  i  OjVJJ 
1NTC35CC 
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FILE:  WVUfi 


F'.RTRAN 


Cl  CuNTIhLE 
02  J» i-i 

COARSE  SCAN  FCk  ZfcRG  PASS 

NEXT  *  LINES  ARE  EXECLiED  IN  FIRST  SEARCH  CNLY 
YC=  <  YC-Y (J  i  >/ ( 1 1 J«i  )-> I J  I » 

ZXIl)*tZil ,o*l 1-- U 
IflZAt  M.EO.HT)  uJ  7C  10/ 

ia*2 

12  JO  1  C3  I*13»NX 

ZXMJMZll.JMI-ZlIrJ) )*  YC+Z I 1 »J)  .  . . 

ZM~M(i  »«i*-iJ-Zli-iiJ  J 1  *  I Y  G-*Y  ( J  ))/( Y< J+ll-YI JI)*Z4I-1»JI 
IF  U«(l  l.cO.rtfl  oC  (C  ICd  „ 

IF  KZXl  I-U-hTi/ IZXI  11-PT)  .LT.o.O)  GO  TO  10* 

03  CJNTINJc 
GO  TO  HI 

CALCJL  AT  10*1  OF  OEOETRATHNS  IN  EITHER  TRIANGLE 
C4  XHO  (1-1  J»(xt  i  )-x  <1-1  J  )*<YG-Y(  JJI/  IY  I  J*1  J-Y  <Ji.J 
IFlI'l.cO.HlI  GJ  'i  3  10  j 

IF  I  UX  m-HT  1/ UV-FT  ».LT  .J.OI  GOTO  103 
DL)-*Af'  S(iU-i.v)+Ab<:  (ZX  U-l)-nTl 
XS=XI1-IHIX<-*I  I-l)  )  *  io  j(  Z  a(  l— D-HT)  /DO 
GC  TC  I  i  0 

C5  >V=>N-MxUJ-X,X»*A3j(ZH  -HT1 /(  AbSUM-HT)+ABSUX  (ll-HT  li 
GO  TC  llo 

LCulCAL  F  ti\  t  T  R  AT  IONS  ACCLPULAT  1CN 
01  XRi  1  J*=x(  iJ 
1*1 

GO  TO  113 

oa  xR<P)=x<n 

GO  T C  111 
U  XR< I R|*  >i 
11  lFIIP.tO.ZJ  GO  TO  114 

13  IR=2 
IB- Ml 

GO  TC  112 

1C  IF  1  lR.EC.D  oC  TC  117 
XR I  2  J* lOOJ.U 
CO  TC  11* 

11  >Rl  11  =  1000.0 

XRI2I-1CUJ.J  .  . 

14  PETUPN 


GO  TO  114 


SUBROUTINE  MTPlTIX,  V.Z.XG.  VG.NX.  IY.ZPES) 

C  L !  i\E  a  P  IMtRFCLATlCN  Ct  CALIF3  AT  SOLACE  TO  EVAL 
C  CP  VALLE  AT  C J-ORO 1  iU Tc S  XG.Yj.  Tn t  RESULT  R  ETUKNtG 
Cl  yERS  j  .1*1  X  MO)  »Y  <-»G)>Z(*0»40) 


LJATE 
G  IS  Z  r  ES 


SEARCH  FCK  l  LF  LEFT  X  LIRE 
DO  1  1C*1,NX 
1  F  ( X  »  1C  J  .CL  .XG  1  u*  TC  2 
CONTI MCE 

l*JC-l  ,  „ 

SEARCH  FCP  J  Cr  LCWEP  Y  LINE 
DO  3  JS=lt  IY 
IF  «Y  IJO  .tC.YGJ  CO  TO  4 
CCMIRLc 
J-JC-1 

AL* (YIJ  +  1J-Y (J 1 M  IX C  1  +  11— XC  I ) ) 
BL*  Y ( J) -A  L*  X  til 
7CR-  I  AL *X0 i ♦3L 
IFlYG.cC.YC3)  uj  TT  I 
IFIYG.il.YCU  G  J  TO  5 
RESULT  If.G  l  IS  II  L  J «R tR  TRIA1GLE 
X3*X  <1+1  ) 

-Y3*  YIJ  1 
2J-Z  l  HUJI 
GO  TO  6 

RASULT  l‘l  C  Z  IS  IN  UPPER  TRIANGLE 
X3»X  U) 

Y3*Y( J* 1) 

23*  Z  1  It  J'«  1 ) 


INY03510 
1NTU-520 
INTC3530 
INT  03S40 
I N  T  o  3  5  S  J 
iNTOiSuO 
INT0357J 
INTC35EO 
IITC35SO 
lNTOjt.00 
INTC341C 
lNTOjtZO 
INTuiodO 
IN  TC3E4G 
INTO  16  SO 
INTG3ct.ll 
INTCit IJ 
I  NT 03aOo 
INTCJ6S0 
1NTo32CG 
INTG37U 
I ITC372C 
INTO  3  730 
IN  103740 
INTCJ75G 
INT  03760 
IN TC377o 
INT  037  80 
I  NT 03  7  40 
IN  TC38CO 
INT036IO 
INTO 3820 
IN  T03d  30 
INT 03c 40 
INT038S J 
INT C3 EtC 
INTtMfc 10 
INTC3tEC 
INT  038  90 
I N  T03900 
INTC3S10 
INVC3S20 
I N  TO  j'J  30 
INTC3S4C 
INI  0  '34  SO 
I N  TO  3StO 
INTC2S7G 
INTO 39 80 
InT  C3SSC 
IN  T  0*000 
INTO* 0*0 
INTC4C20 
IN  T0*0 30 
I  NTJ4040 
IN  TC* J 30 
INT  J*OfcO 
I N  T  oi  0  7  3 
INT  C*CeO 
1 NT  0*0 90 
INTC41C0 
INTO*  MO 
I N  TC* 120 
me*  i3o 
INTO*  1-*j 
INTJhISO 
INTO*  MO 
1  NY  0*1  '0 
1 N  TO* 1  SO 
INT  C* 1  SO 
INT  04200 


77 


fro 


FILE:  WVC3  FORTRAN  Pi  HP  S 

GENERAL  CAlCLLaTI  J.N  :«CjU  FCR  bCTN'T*!  ANGLES 
Xl>X(li 
Vl*ViJi 
2 1*2 ( It  J  ) 

Xi*XUU  I 
T2*Y<  J*U 
li*l ( I ♦ l , J* LI 

PP*lUl-<Ji*lYl-Y*i-tZl-ZZi*lYl-Y3l  >/nXl-X3)*|Yl-Y2)-I  X1-X2) 

iBi« iln-L  J-SP*(  Xi-Xi  U/t  Yl-lul 
CP*ZI-AF'M1-^P*Y  i 
2R  ES»AP*  >(i*rtP*  »i>*CP 
60  TO  1JJ 

C  RESLLTING  l  IS  CN  JIVUIAG  LINE  .  ... 

7  EXES*  IE  I  1»1,0*1J-Z(  1  ,  J  ) )  *1  X^-XII  I )  /(  X(IMI-XI|il»aiijl 
100  RETUhF 
EM  J 

SU^CUT  1NE  M  Ax AT  Y ( £  »  Z'liX  »  IX.NY  )  „  , 

C  EVALlATi  CM  UF  f* AX  JF  CF  CALI  ARAT10N  SURFACE  NATR  IX 

DIMENSION  ZIAO.AL) 

CO  1  1*1. NX 
00  1  J*I .XT 

1  1FIZC  I.  J  I.Jc.EMAX)  rux*m,ji 

RETURN' 

ENJ 

SUBRCUTI'IE  NINaRyIZ.ZMN.NX.NV) 

C  EVALLAT I  C  N  Cr  PiM'M*'  CP  CF  CAL  i  £<  AT  13N  SURFACE  NATRIX 
OIMENSMN  l  IAC.AL) 

00  1  1*1,  IX 
CO  1  J  *  1  t  N  Y 

1  IFIZI  I,  J  I.LE.ZNl.J  »  ZP1N*Z<  l,J) 

RETURN! 

EMU 

C  EMC  OF  WVOT 


INTC<.2  U 
INTo'.ZZO 
1NTCA2  30 
IflTO'KtO 
lMT0»2*.J 
INTL'.^oO 

Iht  jw  /o 
l‘ITCt2eJ 
IN  V  CAZ  SM 
IMCAOGJ 
INT  L7.  J  10 
IN  TO  A  a  20 
1NT0.3J0 
I..TCA2AJ 
1  NTi)-*  j  SO 
1 N  V  L  ■.  3 (• L 

iHTC'titij 
IN  T  C  T  j  SO 
INTOTAOO 
INT  CT410 
iMTJ4-,;w 
l NY JA 4 5 J 
I N  To A A  AO 
INT'C^bO 
INTC'.ASj 
IMTC-.A7C 
INTO<.<.dJ 
INTC'.ASO 
INTCAOCO 
INTO',  ii  10 
INT0*»520 
INT  0>»t»2J 
iNTO'.b'.J 
INT  CAb  50 
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